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1. SOLID STATE 


1. 1 MATERIALS AND SEMICONDUCTORS 

1.1.1 Urate Crystal Nucleation and Growth in Relation to 
Gout and Calculi 

GK- 17042, National Science Foundation 


W. R. Wilcox, A. Khalaf, M. Lin and K. Chen 

Obj ective: 

To gain understanding of factors influencing crystallization of monosodium 
urate, particularly in its relationship to gouty arthritis. 

Approach: 

Nucleation and solubility determinations are being made using a microscope 
hot stage. Impurity incorporation during crystallization is being studied by 
chemical analysis. The crystal structure of uric acid is being determined. 
Solubility and nucleation is being analyzed theoretically. 

Recent Progress: 

The influence of a variety of soluble additives on solubility and 
nucleation of sodium urate have been determined. For example, calcium 
ion was found to greatly enhance nucleation. As pH decreased from 7.3 
to 6. 5 nucleation increased. Nucleation was easier when 5 vol. % of 
synovial fluid from a gout patient was added than when fluid from a 
rheumatoid arthritis patient was added. Cupric ion altered the absorption 
spectra of urate. X-ray diffraction peak intensity diminished as crystalli- 
zation rate increased. The infrared spectra of rapidly crystallized 
"amorphous* 1 urate was identical to slowly crystallized crystalline urate. 

References and Recent Publications 

1. W. R. Wilcox, A. Khalaf, A. Weinberger, I. Kippen and J. R. 
Klinenberg, "The Solubility of Uric Acid and Monosodium Urate", 

J. Med. Biological Eng. (in press). 

2. W. R. Wilcox, "Crystallization Flow", J. Crystal Growth 12, 93 
(1972). 
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3. A. A. Khalaf and W. R. Wilcox, "Nucleation of Sodium Urate in 
Relation to Gouty A r thritis ", 24th ACEMB, Las Vegas, Oct. 31 - 
Nov. 4, 1971. 


1.1.2 Corrosive Effects of Smog 

National Association of Corrosion Engineers 


W. R. Wilcox, T. Shon 


Objective: 

To investigate the influence of smog on metallic corrosion, especially the 
mechanism by which ozone influences corrosion. 

Approach: 


Potentiostatic and ellipsometric measurements are being made during 
exposure of metals to ozone in concentrations typical of smog. 

Recent Progress : 

Having shown that ozone does influence corrosion, an apparatus is now 
being set up to. do potentiostatic measurements in order to investigate the 
mechanism. 

References and Recent Publications : 

1. n The Influence of Ozone on Corrosion", National Meeting of NACE, 
St. Louis (March 1972). 


1.1.3 Ne w Methods for Growth and Characterization of GaAs 
and Mixed UI-V Semiconductor Crystals 

DAHC-1 5-71-G6, Advanced Research Projects Agency (DOD) 

W. R. Wilcox, C. Chang, V. Yip, V. Kuo, A. Esquivel, 

S. Sen, E. Johnson, P. Leung, W. Allred 

Objective: 

To develop new and improved techniques for growth of GaAs and mixed 
III-V crystals in a coordinated program of growth and characterization. 
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Approach: 


GaAs is being grown by a new Czochralski technique, by the travelling 
heater method, and by liquid encapsulated floating zone melting. The 
growth techniques and phenomena are themselves being studied. The 
crystals are characterized by X-ray topography, etching. Hall measure- 
ments, cathodoluminescence, etc. 

Recent Progress : 

GaAs crystals are being routineley grown by the new liquid-seal 
Czochralski technique. Mobilities are too low, so experiments are 
planned with different crucible materials. Up to two passes have been 
made in the liquid encapsulated floating zone melting of GaAs, but 
experimental difficulties are being encountered. A new casting technique 
has been developed for producing feed rods for the travelling heater 
method. Heat transfer in the gradient freeze and travelling heater 
methods has been treated theoretically. Experiments on organic analog 
have enabled evaluation of some of the parameters. 

The maximum rate at which solid particles can be pushed by a 
growing crystal was found to increase with increased stirring. Particle 
mixtures were separated by using a gradually increasing freezing rate. 

Hall measurements on bent GaAs samples provided evidence that 
dislocations may reduce mobilities. Slip bands were evident when bent 
samples were examined by cathodoluminescence in the electron microprobe. 

References and Recent Publications: 


1. R. J. Baughman, R. A. Lefever, W. R. Wilcox, '‘Evaporation of 
Sodium Chloride Melts ”, J. Crys tal Growth 8 , 317 (1971). 

2. W. R. Wilcox, "Movement of Liquid Inclusions by Centrifugation", 

J. Crystal Growth 13/14, 787 (1972). 

3. C. E. Chang and W. R. Wilcox, "Vitreous Boron Oxide: Drying and 
Moisture Absorption", Mat. Res. Bull. 6, 1297 (1971). 

4. W. R. Wilcox and P. J. Shlichta, "A Survey of Movies on Crystal 
Growth", J. Crystal Growth. 15 , 61 (1972). 

5. K. Chen and W. R. Wilcox, "Boiling and Convection During Movement 
of Solvent Inclusions in Crystals ", submitted to Ind. Eng. Chem. 

Fund. 
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1.1.4 Charge Storage in Silicon Nitride and Applications to 
Computer Memories 

F44620-71 -C -0007, The Joint Services Electronics Program 
W. R. Wilcox and W. L. Paterson 


Objective: 

To gain improved understanding of charge storage in silicon nitride films 
so as to be able to tailor-make computer memories utilizing silicon nitride 
films . 

Approach: 

Thermally stimulated currents will be measured in chemically- vapor 
deposited films doped with impurities. 

Recent Progress : 

Undoped films of silicon nitride in thicknesses ranging from 500A to 
3 00 0A have been grown on single crystal silicon substrates employing the 
reaction of silane and ammonia under a variety of preparative conditions. 
Capacitor structures employing a counterelectrode of aluminum or gold 
were fabricated for evaluation. 

For the evaluation of the capacitor structures, apparatus was 
constructed for capacitance vs. voltage and thermally stimulated current 
measurements. The C/V apparatus employed phase sensitive detection, 
allowing small signal levels and ready decomposition of the equivalent 
circuit. The thermally stimulated current apparatus provides linear 
temperature ramps at rates 0. 1 to 10.0 degrees from slightly below room 
temperature to 550 deg C. A high degree of conformity to the linear 
temperature vs. time relationship was obtained by feedback control of the 
rate rather than the temperature. 

Computer studies of TSC spectra and their decomposition by a linear 
programming method were performed. Moderately separated peaks can 
be resolved, but broadened peaks and continua of levels have an ambiguous 
interpretation, with especial error occurring in determination of frequency 
factors . 

Trapping levels deeper than 0. 9 electron volt were found to be 
responsible for the semi-permanent (several year) charge storage effect 
which will be exploited for nonvolatile NDRO memory. These levels 
correspond to TSC peaks occurring at temperatures greater than 3 50 deg C, 
which is consistent with the demonstrated ability of these semiconductor 
memory cells to usefully retain data at temperatures in excess of 
1 50 deg C. 
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1.1.5 Solution Growth of Bulk Crystals of Mixed III-V Compounds 
for Electronics Device Applications 

F44620-71-C-0007, The Joint Services Electronics Program 
DAHC- 1 5-7 1-G6, The Advanced Research Projects 
Agency (DOD) 

W. R. Wilcox and V. Yip 

Objective: 

(a) To prepare bulk crystals of mixed III-V compounds. 

(b) To study physics of bulk solution growth by travelling-heater method. 
Approach: 

A solvent zone is passed through a solid rod by means of a moving heater. 
The crystals are examined for perfection. The interface position and 
shape is determined as a function of operating crystals. It has been 
discovered that GaAs crystals can easily be produced by the gradient- 
freeze technique. 

Recent Progress: 

Progress was made in understanding the physics of the travelling -heater 
method (THM) by heat transfer calculations and experimental determination 
of interface shapes and positions and temperature profiles. While 
exploring means to prepare cylindrical feed rods for the THM it was 
discovered that single crystals of GaAs can easily be prepared by the 
gradient freeze technique. These crystals have moderate dislocation 
contents and electrical properties, so are only suitable as substrate 
materials for epitaxial films and not directly for devices. 
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1.1.6 Epitaxial Growth of III-V Films from Oxygen Pumped 
Solutions 

F44620-71 -C-0007, The Joint Services Electronics Program 
DAHC- 1 5-71-G6, The Advanced Research Projects 
Agency (DOD) 

J. M. Whelan 


Objective: 

To evaluate feasibility of reducing the growth temperatures of epitaxial 
films and increasing their impurities. 

A pproach: 

The stability of oxides on the substrates requires very low oxygen activities 
for their removal by the growth solution. We propose to achieve the 
prerequisite activities at lower temperatures by oxygen pumping using 
solid electrolytes. 

Recent Progress: 

An apparatus was completed for growth of GaAs films from Ga 
solutions with low oxygen fugacities. A stabilized zirconia tube, a solid 
oxygen ion conducting electrolyte, was used as the solution container. 
Oxygen was pumped from the gallium solution electrolytically at rates 
approximately one hundred tirre s faster than could be achieved by chemical 
removal with high purity hydrogen at 800°C. This established the first 
feasibility goal. This coincided with a recent downward estimate made 
elsewhere of the minimum equilibrium oxygen pressures which could be 
established with stabilized zirconia. This was judged to be sufficiently 
significant to warrant a redesign of the apparatus so that oxygen fugacities 
in gallium solutions can be measured which are two to three orders of 
magnitude below those of oxygen saturated solutions. It is anticipated 
that this approach may lead to reliable estimates of the oxygen content 
in GaAs samples. If so, this would be useful for resolving the adverse 
and beneficial effects of oxygen in GaAs devices which have been suggested 
over the past 15 years. 

References and Recent Publications: 


1. Dissertation: Joseph E. Greene, "Glow Discharge Spectroscopy for 
the Analysis of Thin Films", Ph. D. , 1971. 
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1.1.7 Luminescence Due to Native Defects in CdS 


F44620-7 1 -C-0007 , The Joint Services Electronics Program 
M. Gershenzon and G. Turner 


Objective: 

To relate photoluminescent properties of CdS to native defect concentrations 
with a view toward development and optimization of small visible light 
emitters and lasers. 

Approach: 

Calculate defect concentrations for a series of CdS samples annealed 
under controlled conditions and quenched rapidly. Using the samples, 
low temperature photoluminescent and absorption spectra are measured 
and then related to the defect concentrations. 

Recent Progress; 

Samples have been rapidly quenched from both high and low temperature 
and photoluminescence remains bright but is different in character. Angie 
lapping experiments indicate quenching is fast enough. Reannealing 
indicates that samples can be cycled from one condition to another and 
returned reversibly. Photoluminescence shows many lines but absorption 
measurements are possible only on the lines of intermediate wavelength. 
Anneals corresponding to the end points of the series have been completed 
and should represent a 20:1 change in sulphur vacancy concentration. 


1. 1.8 Radiative Recombination in the Alloys Ga, In P 
i 1 -x —x — 

and Ga. A1 P 
1 -x — x — 

F44620-71 -C-0007, The Joint Services Electronics Program 


M. Gershenzon and W. Hawk 


Obj ective: 

To determine the effects of alloy composition on the luminescence 
properties of the III-V alloys Ga. In P and Ga, Al P. 

1 -X X 1 -X X 

Approach: 

Alloy crystals are grown from a solution of excess metal and by liquid 



phase epitaxy. Starting with the known luminescence spectra of GaP, 
shifts and weakening of these spectra with increasing alloy composition 
x (x = 0 is GaP) are used to determine changes in alloy optical properties. 

Recent Progress: 

The donor-acceptor pair lines in Gai_ x In x P alloys, which could be 
resolved only for low In dopings in solution grown platelets and epitaxial 
layers at 4. 2°K, have been observed for slightly higher In levels in 
experiments at 2 5°K where the nitrogen bound exciton becomes ionized 
and cannot produce light, resulting in an enhancement of the donor-acceptor 
luminescence. The epitaxial alloy crystals show resolved donor-acceptor 
spectral lines at In dopings where the energy shift due to In alloying is 
large enough to be detected, although the luminescence efficiency is weak. 
By performing accurate and time consuming experiments we will be able 
to measure the energy shifts of these weak donor-acceptor lines and to 
correlate them with the shifts in the sulfur bound exciton. 

The doubly degenerate sulfur bound exciton in GaP is split into two 
levels in In doped crystals, with a splitting that increases as the amount 
of In in the alloy increases. This effect is attributed to a uniaxial strain 
associated with the incorporation of In into the GaP lattice and could be a 
function of crystal growth conditions. 

Crystals of Gai -X A1 X P have been grown from solution in vitreous 
carbon boats in order to avoid the Si contamination which occurs in growth 
inside quartz containers. The luminescence of these crystals shows only 
the broad far pair band of the GaP donor-acceptor pair spectrum, and the 
shift in this band indicates in excess of 10% Al in the alloy crystals. 

Growth of crystals with a lower Al content is made difficult by the large 
segregation coefficient favoring the incorporation of Al in Gai_ x Al x P. 


1.1.9 Radiative Recombination in Silcon Due to Isoelectronic 
Traps and Deep Levels 

F44620-7 l-C-0007, The Joint Services Electronics Program 
M. Gershenzon and K. Kosai 

Objective: 

To search for new radiative recombination centers in silicon with the 
objective of developing efficient near-infrared light- emitting diodes. 
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Approach: 

(a) Investigate indium -doped float zone grown silicon which has been 
lithium diffused for photoluminescence and absorption due to 
lithium -indium close pairs acting as isoelectronic traps. 

(b) Study the luminescent properties of float zone grown silicon doped 
with impurities such as gold, copper, and zinc, known to produce 
numerous deep levels. Also, study the optical effect of diffused 
lithium on these deep levels, since lithium -impurity close pairs are 
expected to form which may act as isoelectronic traps. 

Recent Progress : 

A ^sample from an In-doped crystal was diffused with Li {[Li] ~ 

3x10 cm”3, [in] ~1 )( 10l6 crn"^} and examined using photoabsorption 

and photoluminescence with the sample immersed in liquid, helium. In 
luminescence, radiation due to decay of excitons bound to both neutral In 
and Li was observed, while in absorption only a line due to creation of 
bound excitons at neutral Li could be identified. Observations showed 
that either the expected Li-In pair had not formed, or, if the Li and Iri 
had paired, the complex does not behave as a trap. 

At sufficiently low temperatures (T ^ 10 K) a series of sharp photo- 
luminescence lines were observed on the low-energy side due to decay of 
excitons bound to neutral Li. We have studied similar lines in B- and P- 
doped Si. These lines do not appear at higher temperatures, whereas the 
associated bound-exciton emission can be observed to T > 3 0 K. These 
lines have been previously observed in B-doped Si and it was suggested that 
they were related to the condensation of an electron-hole metallic-like liquid. 
We are trying to determine the exact mechanism responsible for the 
emission lines. 

Wafers of Si were diffused with Au and Li. In those samples containing 
only A u, all of the intrinsic photoluminescence was quenched and no new 
radiation was observed. Silicon with both Li and Au impurities show an 
entirely new emission band having a sharp peak at 0.76 5 eV and an 
associated low-energy broad band. This emission was not observed in Si 
doped only with Li. We believe that Li-Au pairs are behaving aa iso- 
electronic traps. 
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1. 1. 10 Studies of Semiconductor Surfaces 


F44620-71-C-0007, The Joint Services Electronics Program 


D. B. Wittry, K. Y. Chiu and P. A. Sullivan 


Objective: 

To obtain a better understanding of the phenomena that occur when electrons, 
photons or ions impinge on a solid surface in a gaseous environment. 

Approach: 

Measure surface potentials by electron mirror microscopy and by 
secondary electrons using a scanning electron probe; correlate surface 
potential changes with changes in surface recombination velocity; vary 
the nature of the surface states by changing the ambient gas or by 
implanting ions. 

Recent Progress: 

Constructed and tested an electron spectrometer for use with 
secondary electrons in the electron probe microanalyzer. A pape r has 
been accepted for the 7th National Conference on Electron Probe 
Microanalysis describing this spectrometer. 

The first scanning electron mirror microscope images were obtained 
and the theoretical calculations for the deflection bridge were verified. 

References and Recent Publications: 


1. "Techniques for Semiconductor Research with Electron Beam 
Instruments", Sixth International Conference on X-Ray Optics and 
Microanalysis, Osaka, Japan, Sept. 5-10, 1971. 

2. "Studies of Semiconductors Using Optical Signals from Electron 
Microprobe Instruments ", Symposium on Scanning Electron 
Microscopy, Japan Society of Applied Physics, Tokyo, Japan, 

Sept. 14, 1971. 

3. "An Electron Spectrometer for Quantitative Measurements of Surface 
Voltages with Electron Microprobe Instruments", K. Y. Chiu and 

D. B. Wittry, accepted by the 7th National Conference on Electron 
Probe Microanalysis. 

4. "Applications of Electron Microprobes to Solid State Electronics", 

D. B. Wittry, in Microprobe Analysis, ed. by C. A. Anderson, 

John Wiley Sons, N.Y. (in press). 


- 10 - 




5 . 


"Crystallographic Information from Scanning Electron Microscope 
Images", D. B. Wittry, Proc. Sixth International Conference on 
X-Ray Optics and Microanalysis {in press). 


1.1.11 Electron Optical Characterization of Semiconductor 
Materials 

F44620-7 l-C-0007, The Joint Services Electronics Program 


D. B. Wittry and R. E. Gauldin 


Objective: 

Provide service work for other JSEP research; perform research 
involving new techniques for characterizing semiconductors in the 
transmission electron microscope. 

Approach: 

Use the electron probe microanalyzer for x- ray microanalysis ; use 
selected area electron spectrometry in the transmission electron microscope 
to obtain information about composition and chemical bonding in the TEM 
specimens. . 

Recent Progress : 

Installed and reconditioned the Hitachi HU10 instrument, constructed 
adapters for the magnetic prism spectrometer, made preliminary tests 
of the spectrometer and developed techniques for making thin semiconductor 
specimens. A system for feedback regulation of the high voltage supply of 
the microscope has been designed and is now under construction. 

References and Recent Publications: 

1. "Modification of the High Voltage Supply of an Electron Microscope 
for Improved Stability", R. E. Gauldin and D. B. Wittry, paper 
submitted to the 13th Annual Meeting of the Electron Probe Analysis 
Society of America. 
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1. 1. 12 Cathodoluminescence of GaAs and Its Alloys 

DAHC - 1 5- 71 -Gb, Advanced Research Projects Agency (DOD) 


D. B. Wittry, W. Lin, H. C. Marciniak 


Objective: 

To study mechanisms for spontaneous radiative recombination and 
conditions for stimulated emission of electron beam excited GaAs and 

GaAs x P l-x- 

Approach: 

Observe cathodoluminescence spectra of various temperatures and 
electron beam current densities on GaAs and GaAs-GaP alloys near the 
direct-indirect crossover point; introduce assymmetry by strain and by 
magnetic fields; try to correlate results with theoretical models for 
recombination processes. 

Recent Progress; 

Cathodoluminescence spectra have been obtained from several n type 
specimens of GaAs and GaAs-GaP alloys as a function of current density 
in the electron beam. No evidence for super-radiance has been seen at 
room temperature. Attempts have been made to reduce surface 
recombination by a surface alloying technique, however, these have also 
not been unsuccessful as yet. Additional experiments are in progress on 
specimens subjected to uniaxial compression. 

The results at temperatures down to 30°K on GaAs-GaP alloys are 
encouraging and further experiments as well as theoretical interpretation 
of the results is in progress. 


1. 1. 13 Investigations of Dislocations in GaAs 

DAHC -1 5-71 -G6, Advanced Research Projects Agency (DOD) 

J. Esquivel, D. B. Wittry and W. Lin 


Obj e c ti ve : 

To study the effect of dislocations on recombination processes in GaAs 
after plastic bending at elevated temperatures. 
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Approach: 


Use scanning electron microscopy with detection of infra-red cathodo-^ 
luminescence to observe dislocation densities and configurations; compare 
with x-ray topographs; study both p and n type material with various 
orientations of the bend axis to correlate with electronic properties of 
the dislocations. 

Recent Progress: 

Cathodoluminescence images from n type GaAs (10 /cm Te) and undoped 
(10^/cm^) GaAs have indicated that it will be possible to obtain dislocation 
densities up to about 10^/cm2 with an accuracy of ± 10%. Also, it appears 
to be possible to correlate the dislocation arrays with the known slip planes 
in GaAs. Further studies are in progress using different directions for 
the bend axis and using p as well as n type GaAs. It is anticipated that 
this technique will be an effective non-destructive method of studying 
deformed semiconductor materials. 


1. 1. 14 Avalanche Multiplication in Semiconducto rs 

F44620-71-C-0007, The Joint Services Electronics Program 

C. R. Crowell, C. L. Anderson and Y. Okyto 


Objective: 

To measure ionization coefficients of electrons and holes in Si, GaAs 
and GaP both at low electric fields and in a predicted high field asyptotic 
region which may be useful for stable high gain devices. 

Approach: 

Study of multiplication by photoinjected carriers in Schottky barrier 
configurations --a new approach which should permit a unique separation 
of the effects of multiplication of holes and electrons (very important for 
device (IMPATT) characterization at high frequencies and for noise 
considerations). Theoretical refinement of the way in which ionization 
coefficients are defined is also needed. 

Recent Progress : 

Experimental: 

The development of well-defined Schottky barrier avalanche structures 
is nearly complete. Schottky barrier guard rings in an oxide passivated 
planar structure on lOp. thick, 2 5Q c m epitaxial layers are being used. 
These structures are intended to permit a clear distinction between 
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electron and hole initiated ionization processes. This distinction is 
important for understanding the noise and optimum operation of IMPATT 
and photo-avalanche diodes. 

Theoretical: 

a. Pho toinj ec tion. The above experimental measurements require 
a clear understanding of the electric field dependence of photoinjection in 
Schottky barriers. A theoretical calculation of quantum mechanical 
tunneling effects for this case has been completed for a barrier model 
which for the first time correctly includes effects of conservation of 
transverse momentum and Thomas-Fermi screening in the metal. 

b. Coupled Electron Hole Ionization Rates. Theoretical work 
revealed that the conventionally accepted definitions for electron and hole 
ionization coefficients implicitly convolve effects of electrons and holes 
which depress the coefficients considerably (i.e., according to current 
usage). A modified interpretation allows removing the requirement of 
anomalous threshold energies (too large by a factor of 4) which has been 
needed to explain recent experimental results in the case of Si. 

References and Recent Publications: 

1. ’’Threshold Energies for Elec tron-Hole Pair Production by Impact 
Ionization in Semiconductors", C. L. Andersen and C. R. Crowell, 
Physical Review, B. 5, pp. 2267-2272 (1972). 

2. ’’Energy Conservation Considerations in the Characterization of Impact 
Ionization in Semiconductors", Y. Okuto and C. R. Crowell, submitted 
to Physical Review. 


Experimental Analysis of Impurity Energy Levels in 
Semiconductors 

DA-ARO-D-31- 124-71 -634, U.S. Array Research Office, 
Durham 

C. R. Crowell, C. L. Anderson, M. Beguwala and 
K. Nakano 

Objective: 

To demonstrate how impurity energy levels in semiconductors can be 
characterized by capacitance studies using Schottky barrier (metal- 
semiconductor) contacts. Achieving this aim also clarifies the effect of 
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impurities on the capacitance of semiconductor devices and hence on their 
charge control capabilities. 

Approach; 

The project encompasses both theoretical and experimental work. From 
measurements of the frequency and bias dependence of Schottky barrier 
capacitance and conductance, impurity levels, the capture cross-sections 
for majority carriers, and the concentration of the impurity species will 
be determined. 

Recent Progress : 

a » Theoretical : Crowell and Nakano's theoretical work on the 
frequency dependence of Schottky barrier capacitance is in press. 

Beguwala and Crowell are developing a more comprehensive solution in 
differential equation form which enables one to evaluate the complex 
capacitance (i. e. , conductance and susceptance) of a Schottky barrier in 
the presence of several deep lying impurities. Our approach also facilitates 
the modeling of the impedance with lumped circuit elements. We have 
established that our previously reported three element model is not 
particularly satisfactory for a distribution of deep levels and even for a 
single deep level does not bear a simple relationship to the impurity time 
constant and concentration. Our present approach yields an equivalent 
circuit with 3r+l frequency independent elements for a system with r 
impurity levels. The R^C^n product equals the time constant of the nth 
level over the frequency range of interest. The capacitances have simple 
relationships to the concentrations and energy levels of the different 
impurity species. This model can be uniquely synthesized from experi- 
mental C vs in measurements at a single temperature if one can measure 
capacitance far enough into forward bias. Failing this one need to consider 
C vs temperature and C vs V measurements along with C vs to measure- 
ments for a unique synthesis. A paper on this work is being written for 
publication. 

k* Experimental: Roberts and Crowell have submitted for publication 
their experimental work titled "Capacitance Effects of Gold and Copper 
Impurity Levels in Pt-N Type Si Schottky Diodes". Procedures for 
preparing good Schottky and ohmic contacts to both p and n-type silicon 
are being refined. The electronically balanced broad-band capacitance 
and conductance measurement system is being revised. Pilot measure- 
ments indicated a need for a slower bias sweep to enable use of the longer 
time constant of the phase sensitive amplifier for a higher, signal to noise 
ratio. This modification is under construction. The preamplifier has 
been redesigned and is being rebuilt to operate near critical damping. 

This will allow the system to measure capacitance transients which occur 
within a time span of a few seconds. The system can now display directly 

recorded data in the form of — , and £nC vs V. Facilities for 

c c 3 
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programming C-ti) measurements arc under construction. 

Recent Publications: 

1. c. R. Crowell and K. Nakano, "Deep Level Impurity Effects on the 
Frequency Dependence of Schottky Barrier Capacitance ", Solid State 
Electronics, in press. 

Z. G. I. Roberts and C. R. Crowell, "Capacitance Effects of A u and Cu 
Impurity Levels in Pt-n Type Si Schottky Barriers", submitted to 
Solid-St. Electronics. 


1.1.16 GaAs Growth and Characterization 

DAHC 15-71-G-6, U.S. Army (ARPA) 


C. R. Crowell, C. L. Anderson, S. Joshi, C. Huang 


Objective: 

Task "a": Tunnel and Thermal Effects in Photoemission from Schottky 

Barriers : To determine the corrections to a photothreshold 

plot which is customarily used as primary indication of 
Schottky barrier height. 

Task "b": Schottky Barrier Capacitance Characterization of Impurities : 

To determine the extent to which capacitive characterization 
of Deep Level Impurities is possible in GaAs. 

Task "c": High Impedance Hall System: To construct a high impedance 

temperature controlled Hall system for impurity characteri- 
zation in GaAs. 


Approach: 

Task "a": Theoretical calculations of quantum mechanical transmission 

coefficients of electrons incident in a Schottky barrier are 
convolved. The photoexcited spectrum of electrons are observed 
to predict field dependence of the photocurrent. 

Task "b": Theoretical guidelines provided by work under DA-ARO-D- 

31-124-71-G34 are being used to analyze experimental 
m eas urements . 

Task "c": A balanced varactor bridge amplifier with guarded feedback 

to suppress leakage and common mode effects is being made 
for use in conjunction with a liquid helium Hall effect 
cryostat. 
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Recent Progress: 


Task "a": The paper accepted for publication covers only thermal effects 

in photoemission, but is useful because the normalized 
approach is useful no matter what the metal-semiconductor 
combination. Our extension of the approach with the 
tunneling effects is now near completion and will be useful 
in interpreting results on GaAs where tunneling effects appear 
to predominate over thermal effects at 300°K. 

Task "b": Point-by-point C-V measurements were made on Pt-n 

type GaAs Schottky barriers with frequency dependent 
results observed, particularly near zero bias and in forward 
bias. Further detailed data is required before circuitry is 
completed for C-V and C-W scanning. 


Task M c": The Hall system circuitry is now operational and has 

demonstrated a capability of measuring samples with 
impedances of the order of 10 10 Q with good accuracy. 
Effective leakage impedances of the voltage measuring 
circuitry and the current source are ^10 12 ohms. A low 
impedance mode of operation with microvolt sensitivity 
and digital readout is also available within the same system. 

References and Recent Publications: 


1. "A Convenient Operational Equivalent to the Fowler Photothreshold 
Plot, " by Crowell, Kao, Anderson and Rideout, Surface Science, 
accepted for publication. 


1 - 1 - 17 Experimental Studies of Defect-Induced Vibrational 
Mode Absorption 

F44620-71-C-0067, The Joint Services Electronics Program 
GK-2 5144, National Science Foundation 


W. G. Spitzer 


Objective: 

To accomplish identification and characterization of a variety of point 
defects in semiconductors. In most cases these defects involve impurities 
introduced as dopants in the semiconductors. Once identified, the 
changes in defect concentration as a result of changes in environment 
can be observed directly. 
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Approach: 

The identification is accomplished by measuring the high frequency 
localized vibrational modes associated with the impurity. At present, 
the modes are observed by infrared absorption although it is intended 
to extend the measurements to include Raman scattering. 

Recent Progress : 

a) GaAs: Mg 

The interest in Mg, which is an acceptor in GaAs, results from some 
ambiguities in the interpretation of the infrared absorption of Li- 
compensated material doped with Mn, Cd, or Zn. The models proposed 
for the Li-acceptor ion pair were questioned by some later E.P. R. 
results. The infrared absorption of Mg-doped GaAs compensated by Li 
diffusion is measured and the results are compared with those for 
similar material compensated by electron irradiation and for double- 
doped material. The absorption bands observed in the GaAs: Mg-Li are 
consistent with a model of the following type: 

A) There is only one type of (Mg-Li) pair 

B) Each pair has only one Li paired with each Mg 

C) The Li is in an interstitial position 

D) The (Mg-Li) ir has orthorhombic or lower symmetry 

E) The pair axis is probably in a <r 100 > direction 

b) GaAs Ion Implanted with Si 

The infrared absorption of localized modes for ion implanted Si in 
GaAs have been observed. The implants were at room temperature, at 
energies of 1.0 and 1. 5 MeV, and fluences ^lO 1 ^ ions cm"^. Bands 
previously attributed to Si<"j a and Si^ s defects were observed and the 
effect of isochronal annealing studied. A (SiQ a -Si^ s ) pair band was 
observed after annealing. For most annealing conditions the bands were 
found to either grow or show little change with increasing anneal 
temperature. However, between 500 and 650°C substantial reverse 
annealing was observed for both Si^ a and Si^ g . This reverse annealing 
cannot be the result of pair formation and is being investigated further. 
The electrical data for these high fluences and energies, show no 
apparent correlation with the Si defect concentrations deduced from the 
optical data. The conductivity, even after 900°C annealing, is probably 
controlled by residual damage left in the implanted layer. 

References and Recent Publications : 

1. Disseration, Pat C. Leung, “Infrared Absorption by Free Carriers 
in P-Type Silicon", Ph. D. , 1971. 

2. Dissertation: Lyn H. Skolnik, "Application of Infrared Active 
Localized Vibrational Modes to the Investigation of Impurity 
Incorporation in Semiconductors ", Ph.D. , 1972. 
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3 . Dissertation: Albert E. Cosand, "Infrared Absorption Studies of 
Localized Vibrational Modes and Lattice Modes in Ge x Si 

Ph. D. , 1972. 

4. "Investigation of Te-Doped GaAs Annealing Effects by Optical and 
Channeling-Effect Measurements ", J.A.P. 42 , 3 982 (1971), with 
I. V. Mitchell, J. W. Mayer, J. K. Kung and W. G. Spitzer. 

5. "An Infrared Localized Vibrational Mode Technique for Measuring 
Segregation Coefficients", J. Phys. Chem. Solids 32, 1 (1971), 
with L. H. Skolnik, W. P. Allred and W. G. Spitzer. 

6. "Infrared Absorption of Lattice Modes and the Silicon Local Mode in 
G e xS^l -x -Alloys ", J.A.P. 42, 5230 (1971), with A. E. Cosand and 
W. G. Spitzer. 

7. "Infrared Localized- Vibrational-Mode Absorption of Ion-Implanted 
Aluminum and Phosphorus in Gallium Arsenide", J.A.P. 42, 5223 
(1971), with L. H. Skolnik, W. G. Spitzer, A. Kahan, and 

R. G. Hunsperger. 

8. "Localized Vibrational Modes in Semiconductors: Infrared 
Absorption", Fes tkorperpr obleme XI, Vieweg (Pergamon), 1971, 
pp. 1-44, W. G. Spitzer. 
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1.2 QUANTUM ELECTRONICS 


1.2.1 Frequency Coupling in High Gain Optical Amplifiers 

F44620-71 -C-0007, The Joint Services Electronics Program 


J. H. Parks 


Objective: 

(1) To develop experimental and theoretical techniques for studying 
nonlinear wave propagation in high gain amplifiers, and (2) to relate the 
measured nonlinear parameters to applications in laser communication 
technology. 

Approach: 

To measure the nonlinear coupling between two optical fields travelling 
through a dye laser amplifier at frequencies Wj and w 2 . When these 
frequencies both lie within the dye radiative linewidth, the subsequent 
coupling depends strongly on the difference (w 2 -wi). A tunable dye laser 
will provide the variable frequency fields. The related theoretical calcu- 
lations will proceed from a perturbative treatment of^the atomic interaction 
with two travelling wave fields having wave vectors and h 2> 

Recent Progress : 

The saturation behavior of two optical fields in an absorbing or 
amplifying medium is being investigated. An intense, saturating dye laser 
pulse near 4400 A° and a weak probe field at 4400 A 0 from a c.w. He-Cd 
laser are employed in this experiment. The modulation of the c.w. laser 
absorption in Perylene (liquid phase) as a function of the dye laser wave- 
length allows a direct measurement of the nonlinear coupling properties of 
the absorbing transition. The high power dye laser has been constructed 
with an optical pumping source supplied by a frequency doubled ruby 
laser. This will provide dye laser pulses of ~1 megawatt peak power 
and '—2 0 nanosec widths, tunable near 4400 A°. A theoretical model 
describing the nonlinear coupling between two monochromatic fields 
having wave vectors kj and k 2 has been developed. For the case of a 
Doppler broadened medium, the nonlinear response suggests unique 
interactions between the h#o fields even when the fields travel a relatively 
large angles for which | h x - Tx 2 | ~ A(ju d where Awd is the Doppler linewidth. 

This particular model has not been considered previously and leads to 
interesting results which have applications for large aperature high power 
gas lasers such as C0 2 . 
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References and Recent Publications: 


1. ’’Optical Broadband Amplifier”, J. H. Parks, A. Javan, J. Opt. Soc. 
Am. 6J_, 658 (1971). 

2 . ’’Theory of Noise Amplifiers", J. H. Parks, F. Hopf, J. Opt. Soc. 
Am. 61, 659 (1971). 

3. "Broadband Amplifiers ”, J. H. Parks, Esfahan Symposium on 
Fundamental and Applied Physics, Esfahan, Iran (proceedings to be 
published). 


1.2.2 Experiments with Laser Optical Sources 
GH 32401 , National Science Foundation 


W. L. Faust 


Objective: 

To conduct theoretical and/or experimental studies in three areas: 

1* High Intensity Resonance. In the course of previous work 
involving high intensity illumination of GaP material with 
visible laser frequencies, quenching of visible fluorescence by 
the infrared beam was observed. The current objective is to 
substantiate the assignment of this effect to photoionization of the 
bound-exciton states responsible for the fluorescence and to deter- 
mine the minimum energy necessary for such photoionization. 

2. Excitation Spectra. In the course of the previous work 
described above it was also discovered that fluorescent 
spectra are highly dependent upon the visible frequency 
employed for excitation. Current efforts are now being 
directed to further study these effects and toecploit them 
for assignment of specific spectral features in the 
fluorescence to specific impurities. 

3. Raman Scattering. The objective is to construct a detailed 
theory of the relation between Raman scattering and visible- 
infrared mixing. 
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Recent Progress: 

1. High Intensity Resonance Studies. Fluorescence -quenching has 
been found to occur only if the infrared pulse is properly coincident in 
time with the exciting pulse of visible light. This would seem to eliminate 
simple heating as the mechanism and to indicate a mechanism which will 
be more interesting to study. Further, it was demonstrated that 
fluorescence from impurities previously reported as deeply bound is not 
readily quenched, whereas certain shallow systems are more readily 
quenched. A further objective is to illuminate the material with mono- 
chromatic radiation to determine a minimum energy for quenching. This 
should be the binding energy. 

2. Excitation Spectra Study. We have discovered a prominent case of 
enhancement of fluorescence from a particular impurity in the presence 

of others through resonant excitation. The effect suggests that this is an 
imperative direction for future exploitation of tunable visible dye lasers. 

3. Raman Scattering. This work is at an advanced stage in develop- 
ment of the theory of light scattering from thermal and zero-point infrared 
frequency fluctuations in a crystalline lattice. The models predict with 
some success the intensity versus frequency-including resonance peaks, 
widths, assymmetry, and shifts for varying parameters --both for LO 
excitations and for TO excitations. The LO calculation accommodates 

the contribution of a plasma of charge carriers to the linear dielectric 
response; but, it does not account for scattering from carriers. This 
seems to be adequate for the heavily-damped plasma of GaP. The crystal 
is modeled by Lorentz resonance parameters (determined from linear 
infrared experiments) and by the nonlinear parameters of Faust and Henry 
(measured in an experiment in which coherent visible and infrared beams 
are mixed). 
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1.2.3 Laser Induced Damage in Thin Films and Substrate 


Surfaces 

19628-71-C-0220, Air Force Cambridge Research 
Laboratories 


N. Alyassini, T. Colbert, J. H. Parks 


Objective: 

Structural damage in thin films and at substrate surfaces limits high 
power laser operation and the reliability of antireflection coated optical 
components. This investigation is directed toward an understanding of 
the physical processes leading to material breakdown. 

Approach: 

The development of techniques to measure continuous time evolution 
of damage is a primary approach in these studies. A continuous He-Ne 
laser is reflected from the local region which is undergoing damage 
induced by a high power ruby laser pulse. Studies of the time variation 
in the reflected probe intensity yield the temporal characteristics of the 
damage processes. 

Recent Progress: 

The time evolution of damage to thin films and substrate surfaces was 
observed by monitoring the intensity variation of a He-Ne beam reflected 
from the film or surface during irradiation by a Q-switched ruby laser. 
Monolayered films including ZnS, CaF 2 and NaCl on fused silica substrates 
were studied, as well as uncoated fused silica and LiF substrates. When 
the probe beam is reflected from the exit surface of a substrate sample, 
the reflected intensity was observed to decrease within the timescale range 
of 20-200 nanoseconds to a final value significantly less than that prior to 
damage. When the probe is spatially offset from the damage region, the 
reflected intensity decreases and then returns to approximately the pre- 
damage value within 0-2- 0. 3^1 sec. This reversible effect was observed 
to be associated with a propagating disturbance on the substrate surface 
moving with approximately the speed characteristic of a Rayleigh surface 
wave. The risetime and duration of the plasma which accompanies surface 
damage was also measured. 

The reflected intensity from thin films was observed to increase or 
decrease depending on the film index of refraction. The monolayer film 
samples exhibited reflectance changes on the timescale of the ruby pulse, 
however, an offset probe position indicates this change can occur after a 
time delay of ~ 0. 1 - 0. 2p. sec. 
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1.2,4 Time Resolved Spectroscopy with Pulsed Lasers 


A. Scalabrin, J. H. Parks, R . Brown, L. Singer 


Objective: 

Molecular relaxation processes can be measured in detail by observing 
the time resolved emission spectrum from excited electron states. 

This method is particularly applicable to large organic molecules in which 
many competing decay channels can be separated by combined measure- 
ments of time and wavelength dependence. 

A pproach: 

Intense nanosecond laser pulses will be used to produce a high density 
of molecules in a specific excited electronic state within the timescale 
of the relaxation process under consideration. This provides an 
opportunity to study the spectrum of weakly emitting species and resolve 
the ensuing decay time within the laser pulsewidth. 

Recent Progress : 

A. Prompt and Delayed Fluorescence from Benzophenone 

This research has provided the first complete time resolved spectral 
distribution of the prompt fluorescence and in addition new data supporting 
previously described E-type delayed fluorescence from benzophenone. 

The filtered output from a 50 kW pulsing nitrogen laser was passed into 
1 cm square or 6 mm cylindrical Pyrex cells containing benzophenone 
samples. 

The emission was monitored at right angles to the excitation through 
a cutoff filter (2 M aqueous potassium nitrite, 2 cm path length) to remove 
scattered laser light, and then through a one-half meter Bausch and Lomb 
monochromator onto an RCA IP2 I photomultiplier tube. The signal was 
stored and averaged in a Princeton Applied Research Model 160 Boxcar 
Integrator and (or directly) displayed on a Hewlett-Packard 183 oscillos- 
cope (2 50 MHz). 

The spectral distribution of the fast decay from benzophenone in 
benzene solution at room temperature was obtained with the PAR boxcar 
integrator as follows. The signal from the photomultiplier was repetitively 
fed into the PAR instrument (laser rep. rate, 10-30) where it was stored 
and averaged. A segment of the total signal, at various times after 
triggering of the boxcar , was selected and stored by adjusting the time 
delay T d and sampling gate, T g , of the boxcar . The spectral distribu- 
tion of this segmented signal was determined by the choice of T^ and Tg 
in conjunction with a display of the output from the boxcar on a strip 
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chart recorder as the monochromator wavelengths were scanned by means 
of a motor drive. 

The PAR boxcar integrator used in our experiment can be used with 
Tg as short as 10 nsec or as long as 0. 5 sec while the can be made 
very close to 0 sec. Thus, nsec decays can be sampled and time resolved 
from longer decays. By comparison, mechanical devices, such as the 
phosphoroscope, can be employed only for sampling relatively long-lived 
( >msec) emissions. 

A series of spectra at room temperature at different delay times 
reveals the delayed fluorescence component in the emission from benzo- 
phenone. At = 0, the observed spectrum is a composite of the 
fluorescence and phosphorescence. At > 10 nsec, the prompt fluores- 
cence has ended and the spectrum might be expected to reveal only the 
phosphorescence spectral distribution. However, residual intensity is 
observed as a shoulder at 3800 X. This intensity decays at the same 
rate as the phosphorescence as indicated by a constant ratio for {intensity 
at 3800 j? )/ (intensity at 5250 A). This result is consistent with an E-type 
delayed fluorescence mechanism as the source of the residual high energy 
emission. 

In summary, we conclude that the emission from benzophenone at 
room temperature consists of prompt fluorescence, delayed fluorescence, 
and phosphorescence based on our observations of (a) a very short-lived 
emission from benzophenone ( < 6 nsec) at 300°K and 77°K. (a) a longer- 
lived emission {(isec) at 300°K which is a composite of the phosphorescence 
emission and an emission having a spectral distribution characteristic of 
the prompt emission, (c) a long-lived (msec) emission at 77°K which shows 
only the phosphorescence. 


B. Vibrational Relaxation in the Excited States of Large Molecules 

The large number of vibrational modes in multi-atomic molecules does 
not allow a vibrational level to be easily resolved for the purpose of 
measuring vibrational relaxation processes. The use of short (5 nanosec), 
intense pulses provided by a tunable dye laser offers a unique method of 
measuring these relaxation phenomena. As the laser wavelength is tuned 
to excite a molecule into different higher electronic -vibrational levels, 
the spectral distribution of the ensuing emission becomes very sensitive 
to the gas pressure. In particular, the emitted intensity versus wave- 
length distribution will be independent of pumping wavelength and gas 
pressure only when the excited state molecules have relaxed into the 
lowest vibrational level. At lower pressures {.01-1 mm Hg), the 
emission will occur before collisional relaxation and then exhibit a 
spectral distribution characterizing the higher levels. 

Research underway includes the measurement of vibrational relaxation 
within the first excited state of biacetyl using this technique. These studies 
will be extended to molecules of similar structure, glyoxyl and benzil. 
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1.2. 5 Study of Laser -Irradiated Thin Films 

FI 9628-7 l-C-0220, Air Force Cambridge Research 
Laboratories 


L. G. DeShazer, E. A. Maunders, B. E. Newnam 
and H. R . Owen 


Objective : 

The damage resistance of state-of-the-art dielectric thin-film reflectors 
and antireflection coatings to laser radiation is being studied. Particular 
attention is given to the measurement of a threshold for laser damage to 
thin films . 

Approach: 

A TEM o0 Q-switched ruby laser is used to damage mono-, bi- and multi- 
layered coatings of the materials Ti02» SiC> 2 » ZrO^* MgF 2 and ZnS on 
substrates of glass, rocksalt and spinel. Damage threshold energy and 
power densities are measured with attention paid to dependences on laser 
beam spot-size, film material and thickness, substrate condition and 
incipient scattering of the film. 

Recent Progress : 

The damage experiments were continued during the reporting period 
in studying the variation of the damage threshold energy and power 
densities with adjustments in the film, substrate and laser parameters. 
Emphasis was placed on the establishment of a sensitive and reliable 
criterion for theonset of damage. Many investigators of laser-induced 
damage to surfaces of transparent dielectrics have observed a luminous 
surface spark whenever physical damage occurred. Some researchers have 
proposed that the spark is responsible for the observed damage, and it has 
thus become a convenient practice to identify the spark threshold as the 
damage threshold. We have used three methods to monitor the onset of 
damage: a) spark detection using a sensitive photomultiplier with a narrow 
band interference filter, b) laser-induced increase in the scattering of a 
He-Ne laser probe, and c) optical microscopic observation of film 
breakup. 

The interesting result is, that for monolayer films, an increase in 
He-Ne light scatter was observed prior to or at the spark threshold, 
depending on the film material. If the initial He-Ne scatter levels were 
low enough, we would expect the same result for multilayer films also. 

It was shown that lase r-induced scatter occurs before spark initiation for 
films with high indices of refraction. 
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Another recent aspect of this study has been the investigation of a 
correlation between scattering {before damage) from thin-film coatings 
and their damage thresholds. The scattering was measured by the 
intensity of light scattering at about 13 5° from the direction of the 
incident beam. Since there is a strong wavelength dependence for the 
diffuse reflectance from a dielectric film, we used a severely attenuated 
output from the single-mode Q-switched ruby laser as the measuring 
beam. 

References and Recent Publications: 


1. B. E. Newnam and L. G. DeShazer, "Spot-Size Dependence of 
Laser-Irradiated Damage Threshold of Dielectric Thin Films, " 

IEEE JQE 8, 611 (1972). 

2. L. G. DeShazer and B. E. Newnam, International Conference on 
Thin Films, Venice, Italy, May 15-19, 1972. 

3. L. G. DeShazer and B. E. Newnam, "Study of Laser -Irradiated Thin 
Films, " 1972 ASTM-NBS Laser Damage Symposium, Boulder, 
Colorado, June 14-15, 1972. 


1*2.6 Transient Behavior of Transverse Modes in High-Power 
Pulsed Lasers 

AFOSR -7 1-2066, Air Force Office of Scientific Research 
L. G. DeShazer, B. E. Newnam and G. L. McAllis ter* 


Objective: 

We are investigating theoretically and experimentally the transient 
behavior of transverse modes in high-power Q-switched laser oscillators. 
This research will cover an important deficiency in the present charac- 
terization of "single-mode" pulsed lasers emitting high powers, a 
characterization that is urgently needed in order to achieve the total 
capability of giant-pulse lasers in optical systems. 

Approach; 

During this program, a computer analysis will derive for high-power 
Q-switched lasers the transverse patterns, which will be compared to 
experimental patterns of a giant-pulse ruby laser. Emphasis will be 


Northrop Corporate Laboratories, Hawthrone, California. 
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placed on understanding the single mode performance when various 
method of Q-switching (passively or elec tro-optically) are used to pulse 
the laser. 

Recent Progress: 

The computer program used prior to this grant was rewritten and 
upgraded during this report period. The number of radial increments 
was increased from eleven to forty-one, while the number of temporal 
increments per round trip was kept the same as before. This improve- 
ment enabled us to calculate the transverse patterns for cavity Fresnel 
numbers up to ten with single pass gains of twenty. The computer 
program calculates power density and relative phase as a function of 
radial position and time in the aperture plane, the total power (spatially 
integrated) versus time, and the radial distribution of the energy density 
in the aperture plane. Further, these quantities can be calculated at any 
distance from the aperture. 

At present, we have observed that the pulsewidths are the same up to 
Fresnel Number of 5. At Fresnel Number 8 the gain saturation perturbs 
the profile greatly, producing spikes in the temporal profile. Such 
temporal spikes have been experimentally observed with this calculated 
behavior, but this is the first calculation of this observation. 

Several attempts were made during the computer analysis to include 
the effect of a saturable absorber on the transverse mode distortion, but 
several inconsistencies in the calculations resulted. These inconsistencies 
were traced to the method used for the interaction of the forward and 
backward traveling portions of the light pulse with the absorber. Now 
we are exploring several suggestions as to a method of including the 
saturable absorber in the computer analysis. 

References and Recent Publications: 


1. B. E. Newnam and L. G. DeShazer, "Criterion for Obtaining 

Single-mode TEM QO Operation of High-Powered Pulsed Lasers, " 
to be presented at the Meeting of the Optical Society of America, 
San Francisco. 
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1.2. 7 Advanced Nonlinear Absorption Devices for Laser 
Sys terns 

NGL -0 5-018- 044, National Aeronautics and Space 
Adminis tration 


L. G. DeShazer and J. W. Austin 


Objective: 

We are investigating theoretically and experimentally the behavior of 
optical resonators containing saturable absorbers. Even though saturable 
absorbers have been used in laser systems as passive switches and 
modulators since 1964, the operation of these absorbers in conjunction 
with optical resonators has not been adequately described as yet. 

Recently, a comprehensive theoretical treatment of their behavior was 
developed and now this general theoretical approach will enable us to 
design a second generation of saturable-absorber optical devices. 

Approach: 

The features of this problem are succinctly contained in the problem of 
a Fabry-Perot resonator containing a saturable absorber. Such a device 
is called a saturable resonator. A general steady-state solution to the 
performance of the saturable resonator has been obtained and was 
published in the dissertation of J. W. Austin. Now, this solution has 
given suggestions as to a method of including the saturable absorber in 
a computer analysis of the transient solution. 

Recent Progress: 

Presently, we are investigating the time -dependent analysis of a 
saturable resonator. An analytical solution of the rate solutions is 
tractable only for a step- or square-pulse input; hence, for more 
complicated pulse shapes, such as a Gaussian shape, we resort to 
numerical techniques. We have begun an analysis of the saturable 
resonator for pulses having widths much larger than the response time 
of the absorber. Transfer curves describing the transmission of a 
laser pulse have been obtained for this case. For a given peak intensity 
and a specific absorber system the intensity level at which the saturable 
resonator turns "on" can be adjusted by varying mirror reflectivity and/or 
low level absorber transmission. 

By these adjustments significant narrowing of a pulse can be achieved. 
Also it was demonstrated that the transmitted pulse shows an asymmetry 
about the peak; this asymmetry is due to the hystersis of the saturable 
resonator when it operates in the bistable region. Experiments are being 
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conducted to verify these theoretical predictions. 

In addition to these pulse studies, an analytic solution is being 
attempted for the case of a Fabry-Perot resonator with an absorber 
whose transmission decreases with increasing intensity. This solution 
may be of interest for studying laser damage to thin films. 


1.2.8 Laser and Spectral Properties of Rare Earth Crystals 

AF-AFOSR -69- 1622A , Joint Services Electronics Program 
Gift, Union Carbide Corporation 


L. G. DeShazer, J. K. Guha and L. M. Hobrock 


Objective: 

We are developing methods to allow comprehensive application of the 
theories of relaxation processes in rare-earth-ion/solid systems to 
quantum electronic devices of current interest. It is well known that 
certain rare-earth crystals, such as Nd/YAG, have achieved relatively 
high performance as solid state lasers. However, development of rare- 
earth lasers during the past ten years has proceeded primarily on an 
empirical basis. Now, we recognize that the theories of RE 3+ /solid 
systems are sufficiently precise to provide quantitative guidance in 
selecting new, potentially useful RE^+/ solids for use in laser devices. 

Approach: 

We will develop computer software to permit efficient computation of 
radiative transition matrices for all 4f n systems. Concurrently we will 
experimentally determine multiphonon parameters for important host 
materials by performing lifetime measurements on selected RE^+/solid 
systems, and determine the intensity parameters by measuring intensites 
of these selected systems. 

Recent Progr es s : 

There has been a recent shift of interest from YAG to yttrium 
orthoalumina te (YAP) as a host for neodymium. The important difference 
between the two materials is that while YAG is cubic and optically iso- 
tropic, YAP is orthorhombic and optically biaxial. The anisotropic 
nature of YAP makes it of interest because the wavelength, polarization, 
and output energy is a function of crystallographic orientation. This 
feature has been emphasized strongly because it allows certain laser 
characteristics to be optimized through the proper choice of rod orienta- 
tion. We have studied the polarization properties of the NdrYAP spectra 
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as a function of crystal direction. The polarization of each transition 
has a different elliptical nature, but the anisotropic optical properties of 
YAP have confused the analysis of the polarization measurements since 
the birefringence of YAP modifies the polarization states of the spectral 
lines. Our polarization analysis will be completed soon. 

An extensive study of the Nd:YAP laser was made. The input energy 
versus output energy was measured for the two strongest transitions 
(1. 0645 and 1. 0795|a) for various crystal directions. The experiment 
involved a plane mirror resonator with a long pump pulse (600psec). 

The slope efficiency was found to be highest for the a-axis rod with 
chromium, 3.3 times that of YAG. The a-axis rod at 1.08g. with 
chromium seems to be the best choice as a laser material for pulsed 
laser action. The same is not true for CW operation; our preliminary 
studies indicate that the b-axis rod with chromium seems to be the most 
suitable. 

References and Recent Publications ; 

1. L. M. Hobrock, L. G. DeShazer, W. F. Krupke, G. A. Keig and 
D. E. Witter, "Four-Level Operation of Tm:Cr:YA10, at2.35p, n 
IEEE JQE 8, 533 (1972). 

2. J. K. Guha, L. G. DeShazer, G. A. Keig and D, E. Witter, 
"Polarization and Crystal Direction Effects on Laser Performance 
of Neodymium/ Yttrium Or thoaluminate, " IEEE JQE 8 , 53 6 (1972). 


1<2 *9 Theory of Electromagnetic Scattering Processes at 
Resonance, and With Intense Fields 

F44620-7 1-C -0007, The Joint Services Electronics Program 
R. W. Hellwarth and J. H. Marburger 

Objective; 

To be able to interpret and predict the spectrum of the electromagnetic 
radiation that is scattered by matter from an incident, monochromatic, 
electromagnetic wave as its frequency is tuned through a resonance of 
the scattering medium and to understand the modifications of the spectrum 
when the incident field strength is very large. 
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Approach: 


Our approach will be theoretical but directed toward meeting the needs of 
the experimenters who will soon be studying the details of resonance 
scattering for the first time with new tunable laser sources. (Some of 
these sources are also capable of producing very high optical intensities.) 

Recent Progress: 


This is a new task. The work has only recently become of interest as 
practical and reliable coherent infrared and optical sources have 
approached development. 


1.2.10 Investigating Mode Competition Effects in High Power 
Gas Lasers 

DA HOI-7 1 -C- 1476, U.S. Army 


W. H. Louisell 


Objective: 

To study resonators with diverging central portion and converging rims 
to spread out mode volume and simultaneously reduce diffraction loss 
for single mode, high power operation. 

Approach: 

The integral equations describing wave propagation in unstable rimmed 
mirrors were solved by computer. Design parameters were found for a 
fundamental mode with large mode volume, low diffraction loss and 
discrimination against higher order modes. 

Recent Progress : 

Theoretical work has terminated and a paper has been written for 
publication (see below). 

References and Recent Publications: 


1. "Transverse Mode Suppression Using Rimmed Unstable Resonators", 
M. Lax, C. E. Greninger, W. H. Louisell and W. B. McKnight, to 
be published. 
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1.2. 11 Studies of Coherent Optical Devices 

F44620-71-C-0007, The Joint Services Electronics Program 


W. H. Steier 


Objective: 

The invention, evaluation, and application of devices for the control and 
detection of coherent electromagnetic radiation in the infra-red and 
optical spectrum. 

Approach: 

a) Study methods to stabilize and increase the tunability of parametric 
oscillator; 

b) Apply the tunable IR parametric oscillator to monitoring and trans- 
mission measurements in the atmosphere. 

Recent Progress: 

A frequency doubled Nd-YAG laser has been built and used to pump 
a LiNbC >3 parametric oscillator. The pump uses internal doubling with 
LilO^ and a novel prism coupler to improve the conversion efficiency. 
The pump produces more than 1000 W in 200 n sec pulses at 5300 % in 
a very stable TEMoo mode. An experimental and theoretical study of 
two pass harmonic generation using the prism coupler and the Nd-YAG 
laser is complete. 

The LiNb 03 parametric oscillator is single resonant and can be 
temperature tuned from 2. 0 to 2. 5|j. Linewidth and stability measure- 
ments are in progress. 

A pulsed ion laser has been built and operated as a possible pumping 
source. Using xenon the laser can produce up to 2000 W pulses on 
several lines in the visible. The output has several transverse modes 
and the laser cannot be mode controlled without significant loss of 
power. Measurements of small signal gain on several lines were made. 
The pulsed ion laser does not appear useful for parametric oscillator 
pumping. 

A theoretical study of angle tuned parametric oscillator pumping is 
complete and shows the limits of angle tuning as determined by the 
threshold increase due to the non-colinear interaction. The optimum 
cavity design for angle tuned oscillators is described. 

References and Recent Publications : 

1. "Threshold Calculations for Angle Tuned Parametric Oscillators", 
R. Basu and W. H. Steier, presented at the Optical Society of 
America Meeting, - Tucson, Ariz. , April 1971. 
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2 . 


"Bandwidth and Threshold Calculations for Angle Tuned Parametric 
Oscillators", R. Basu and W. H. Steier, IEEE Journal of Quantum 
Electronics (accepted for publication). 

3. "Experimental Study of Two Pass Internal Second Harmonic 

Generation Using a Prism Coupling Technique", D. Gonzalez, 

S. Nieh, W. H. Steier (in preparation). 
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1.3 MAGNETISM 


1*3.1 Magnetos trictive Behavior of the Si-Doped 
YIG 

F44620- 71 -C-0007, The Joint Services Electronics Program 
J. Smit and D. R. Mack 

Objective: 

The technical properties of ferrites, such as permeability, coercive force, 
switching time, resonance, microwave losses, etc., depend upon the 
coupling of the magnetic moment to the lattice. The objective of this task 
is to investigate the magnetos trictive behavior of Si-doped yttrium iron 
garnet (YIG). 

A pproach: 

Strain guage measurements are made between 4°K and room temperature 
on a single crystal of Si-doped YIG in magnetic fields applied in various 
directions. The results are analyzed in terms of non-quenched orbital 
moments of the ferrous ions. 

Recent Progress: 

It was known that Si substitutions have a large influence on the magneto- 
striction of yttrium-iron garnet (YIG) at room temperature. We found 
that its effect is at least 10 times as strong at low temperatures, with a 
maximum, especially for the [ 110 ] direction, at about 80°K. Moreover, 
below 50°K a magnetic anneal effect occurs. This could all be quantita- 
tively interpreted in terms of the divalent iron ions introduced by the 
quadruvalent Si ions. These Fe^ + ions have a non-quenched orbital 
angular momentum, which is strongly coupled to the lattice. At low T, 
these Fe^ + ions stay in the neighborhood of the Si^ + ions, but break away 
at higher T. This increases their magnetos trie tion, because the 
anisotropic perturbing fields of the Si ions then become less important. 
Good quantitative agreement between experimental and theoretical data 
for the anisotropy and the temperature dependence of the magnetostriction 
was obtained. Detailed results will be available in references cited below. 

References and Recent Publications; 


1. In preparation: "Magnetostriction in Silicon-Doped Yttrium Iron 

Garnet". 

2. Dissertation: Daniel R. Mack, "Magnetostriction in Silicon- and 

Calcium-Doped Yttrium Garnet", Ph.D., 1972. 
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1.3.2 The Measurement and Analysis of Hyperfine Fields in 
Magnetic Alloys and Compounds 

F44620-71-C-0007, The Joint Services Electronics Program 


J. Smit 


Qbj ec tive: 

To measure internal magnetic moments in solids. 

A pproach: 

Perform nuclear magnetic resonance experiments and measurements of 
the Mossbauer effect for information about the local magnetization and 
chemical composition of mixed ferrites as used in computers and micro- 
wave devices, and of segregated permanent magnet alloys such as Alnico. 

Recent Progress; 

The NMR study on Ni-based alloys was completed and published 
(see References below). The results provided information about the 
distribution of the excess electrons of the solute atoms over the neighboring 
atoms. The distribution appears to be very important not only for magnetic 
properties but also for many other properties such as the mechanical yield 
strength of similar alloys. 

The experimental techniques and types of analysis used in this work 
will be applied to other alloys in the continuation of this project. 

References and Recent Publications: 


1. "Mossbauer Effect in Li-Zn Ferrites", J. Applied Phys. , 42, 2344-48, 
1971. 

2. "Nculear Magnetic Resonance in Ni-Cu, -Zn and -Ga Alloys", Solid 
State Communications, V 9, pp. 1880-1892 (1972). 
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1.3.3 Study of Paramagnetic Impurities in An ti-ferromagnetic 
Crystals 


F44620-71 -C -0007 , The Joint Services Electronics Program 
J. P. Hurrell 


Obje ctive: 

To understand the fundamental microscopic and macroscopic properties 
of materials exhibiting second order phase transitions. 

A pp roach; 

Using the resonance properties of paramagnetic impurities to probe local 
fluctuations, crystal distortions and/or spin canting. 

Recent Progress : 

The well-localized resonance of a paramagnetic impurity (E r ) in 
specially grown antiferromagnetic MnF^* Associated with this impurity 
is a fairly stable charge compensation which must be dispersed before 
further detailed investigations can be performed. Paramagnetic impurities 
have also been observed in TmAsO^ and TmVO^ in their low temperature 
phases where they behave as Van Vleck paramagnets. In the arsenate, we 
have observed a rapid temperature dependence to the linewidth of an 
erbium resonance in the vicinity of the transition temperature which we 
hope to be able to correlate with fluctuations in the TmAsC >4 system: The 
natural defect center observed in the vanadate does not exhibit a similar 
line broadening. Consequently, we feel that the broadening reflects the 
effects of lattice distortions and not magnetic interactions. 

References and Recent Publications: 

1. 1 'Strain -Induced Effects in the EPR Spectrum of Dy 3+ : YVO 4 ", Bull. 
Am. Phys. Soc. 17, 49 (1972). 

2. "Intensities of the Vibronic Sidebands of Sm^ Fluorescence in 

BaF 2 : Theory and Experiment", Bull. Am. Phys. Soc. 17, 130 (1972). 

3. "Nuclear Induced Electron-Spin-Echo Decay in Solids", Solid State 
Comm. 9> 461 (1971). 

4. "Theoretical Survey of the Sidebands of Sm 2+ Fluorescence in BaF 2 ", 
accepted for publication in Physical Review. 
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1.4 METALS, ALLOYS AND LOW TEMPERATURE PHYSICS 

1.4.1 A n Electron Microscope Investigation, of Grain Boundary 
Phenomena During High-Temperature Creep 
F44620- 7 1 -C-0007, The Joint Services Electronics Program 


T. G. Langdon 


Objective: 

The occurrence of grain boundary sliding in polyc rystalline materials 
during high- temperature creep leads to crack formation and premature 
fracture. This research program will investigate grain boundary 
phenomena using transmission electron microscopy. 

A pproach: 

High-temperature constant stress creep tests will be terminated at 
selected strains and the specimens sectioned for examination by 
transmission electron mic roscopy. Emphasis will be placed on the 
formation of substructure and the development of grain boundary dislocations. 

Recent Progress : 

This is a new project, however, the following recent activities have 
occurred in preparation for formal initiation: 

a) Equipment: A 20,000 pound Instron testing machine has been 

obtained which is one of several important items of equipment 
needed for the research. The machine will allow tests at constant 
strain rate in addition to the constants stress to include crack 
formation under dynamic high strain -rate conditions. 

b) Preliminary theoretical work concerning: 

1) The nature of creep mechanisms in hexagonal materials and 
the shape of the strain: time curves. 

2) The significance of grain boundary sliding under conditions of 
low-stress creep, and in particular, the application of 
computer techniques to determine the stress dependence of 
sliding. 
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1.4.2 Experimental Investigation of the Electronic Structure 
of Intermetallic Compounds and Transit ion Metals 
F44620-71 -C-0007, The Joint Services Electronics Program 

M. H. Halloran, J. H. Brewer and R. D. Parker 


Qbj ective; 

To carry out experiments on: A uln 2 to determine the similarities and 
differences of its electronic structure as compared to AuGa 2 ; on MgZn 2 
to determine its electronic structure; and, on vanadium to determine the 
extent of its electronic similarity to Nb and Ta. 

Approach: 

The deHaas-van Alphen effect will be used in a superconducting solenoid 
at temperatures as low as 1. 1 K and magnetic fields as large as 100 kOe. 
Using uniaxial stress up to 2 kilobars when appropriate, the Fermi 
Surface topologies will be determined for the materials. 


Recent Progress; 

The experimental research on AuGa 2 and Auln 2 is essentially 
complete, and one publication appeared during the last year. The Fermi 
surfaces of A uln 2 and AuGa 2 are very similar with a few very significant 
differences which are consistent with the modifications to the nearly Free 
Electron Model which have been calculated in an APW band structure 
calculation. All the observed sections of Fermi surface can now be 
attributed plausibly to the calculated models, which in turn yield a possible 
explanation for the anomalous temperature dependence of the magnetic 
properties of AuGa . 

w 

The experimental results on vanadium are in progress and should be 
completed in mid 1972. Early results confirm the expectations that 
vanadium should have a Fermi surface similar to niobium and tantalum. 
Experiments have been performed on dendritic crystals grown by the 
Bureau of Mines and on zone-refined crystals grown in our laboratories, 
with the latter giving slightly better results so far. 


The results on MgZn 2 are very preliminary at this time. Samples 
have been grown with resistance ratios of approximately 200 

( ^ )» and very large signals have been obtained using the dHvA effect. 

vk 4.2 
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References and Recent Publications: 


1. ’’The Fermi Surface of Vanadium", by R. D. Parker, J. C. Abele 
and M. H. Halloran, Bull. Am. Phys. Soc. 16, 838, August 1971, 
Seattle, Washington. 

2. "New Fermi Surface Observations inAuIn 2 "> by J. H. Brewer and 
M. H. Halloran, Bull. Am. Phys. Soc. 17, 92, January 1972, 

San F rancisco. 

3. Articles Published: 

a) "New Fermi Surface Results in AuGa 2 M , by J. C. Abele, 

J. H. Brewer and M. H. Halloran, Solid State Comm. 9 , 977 
(1971). 


1.4.3 Thermoelectric Properties of New Materials in 
Metastable Crystalline and Amorphous States 
F44620-71-C-0007, The Joint Services Electronics Program 

R. Wang 

Objective: 

The synthesis and characterization of new amorphous and metastable 
alloys with high potential of good thermoelectric performance. 

Approach: 

To synthesize new amorphous solids and metastable crystalline phases 
by rapid quenching from liquid state; to measure the electronic and 
thermoelectric properties of s true ture- sensitive properties. 

Recent Progress: 

A number of tellurium -based semiconducting amorphous materials 
(Te-In, Te-Ga, Te-Ge) have been prepared by rapid quenching alloys from 
the liquid state. All amorphous phases are characterized by x-ray 
diffraction and electron microscopy. 

The thermoelectric power of an amorphous In 20 Te 80 allo Y was 
measured from 80°+ to 300°K. At high temperatures electrons are excited 
thermally across the gap and the Fermi level will lie near the middle of 
the gap. The thermoelectric power measurements show the intrinsic type 
increasing as temperatures decrease. 
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After annealing at 100°C for 16 hours the samples become crystalline. 
The thermoelectric power then decreases an order of magnitude, and the 
high temperature intrinsic region shifts to higher temperatures. In some 
cases the low temperature drag effect occurs. 

Electric resistivity and Hall mobility measurements are in progress 
but accurate data on the Hall measurements will be difficult to obtain due 
to amorphous low mobility carriers. 

We have found, for the first time, a memory switching effect in the 
In-Te, Ga-Te and Ge-Te amorphous semiconducting materials. The 
switching characteristics and mechanism are now under investigation. 

References and Recent Publications : 

1. In preparation: "Thermoelectric Properties of Te-Based Amorphous 
Semiconductors ", 


1.4.4 Properties of Electronic Materials and Devices a t 
Low Temperatures 

F44620-7 l-C-0007, The Joint Services Electronics Program 
T. I. Smith, Y. B. Kim and M. D. Daybell 

Objective : 

To study the electronic properties of pa rticula r materials at temperatures 
of 4°K and below, with the intent of utilizing these materials in instrumen- 
tation and communications devices operating at these temperatures. It is 
expected that significant improvement in many such devices will be 
realized by operating at low temperatures. 

Approach: 

The measurements fall into two categories. In the first, data will be 
obtained through measurement of the microwave surface impedance of the 
materials, using cavities and high frequency transformer configurations. 

In the second, we will measure the properties of conventional solid state 
devices operated at low temperatures. Quantities such as device noise, 
gain, and operating power level are important. 

Recent Progress : 

The initial efforts in this new project involve measuring the loss tangent 
of very low loss dielectric materials needed to support devices inside 
high Q, superconduction microwave cavities. 
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1.4.5 The Effects of Electric and Magnetic Fields and Other 
Phy sical Parameters on the Nucleation, Structure, and 
Residual Properties of Vapor-Deposited Metal Films 
AT (04-3)- 1 13 , Project 22, Atomic Energy Commission 

L. E. Murr, O. T. Inal, M. C. Lin, H. P. Singh, 

H. R. Vydyanath, and G. I. Wong. 

Objective: 

To determine whether electric or magnetic fields applied parallel to or 
perpendicular to a substrate during vacuum -deposition can significantly 
alter the nucleation and growth of their metal films. 

Recent Progress: 

This 3 year project was completed at the close of the reporting 
period. Detailed results are contained in the cited publiscations 
summarized as follows. 

During the past year, the effects of electric fields on the nulceation 
and growth of In, Sn, Bi, Pb, Cd, Au, and Ag thin films was continued, 
with particular emphasis on the low-melting point metals. No effects 
were observed with carefully defined potential gradients ranging up to 
3 KV/cm in the plane of the substrate. During the course of these 
investigations, several prominent nucleation characteristics were 
observed. These included vapor -liquid- solid (liquid nuclei supercooling) 
and vapor-solid nucleation characteristics which were observed to be a 
function of substrate temperature and film thickness evaporated at 10" 
Torr vacuum. Detailed analyses were made of the effects of substrate 
temperature, evaporation rate, and film thickness on nucleation, film 
growth, and continuity for A u, Ag, Pd, Pt, In, Sn, and Al. The results 
indicated the substrate temperature /melting temperature to be a critical 
factor in film growth, and to influence film continuity as a function of 
mean foil thickness, in addition to the surface energy of the nuclei, and 
the energy of the substrate/overgrowth interface. 

A special apparatus was constructed for the application of a magnetic 
field in the plane of a substrate, having a linear field capability of 0-4 
Kgauss. Films of Ni, Co, and Fe grown under these field conditions at 
a range of substrate temperatures showed no effect of a magnetic field 
in a vacuum of 10"^ Torr. 

Detailed examinations of Pt overgrowths on Mo W Ir field emitter 
end form substrates by field-ion microscopy indicated the very early 
stages of nucleation and film growth to be characterized by an accommo- 
dation of the overgrowth specie to the crystallography of the substrate 
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(basal-plane pseudomorphism). Evidence obtained suggests that the 
overgrowth reverts to its habit crystallography at some critical thickness 
(usually only about 10-20 atom layers in mean thickness), and this 
relaxation can be accompanied by the introduction of crystal defects. 

This inves tigation has shown that the enhancement of thin film growth-- 
viz. , grain size control, the production of large grain or single-crystal 
thin films, is not significantly influenced by an electric field applied 
parallel or perpendicular to the substrate during evaporation for linear 
and reasonably homogeneous field strengths up to 10^ volts/cm. The 
dominant parameter influencing film growth is observed to be the substrate 
temperature, particularly the substrate temperature /melting temperature 
ratio. No significant enhancement of grain size in ferromagnetic thin films 
has been observed for field strengths up to 7 Kgauss applied perpendicular 
to the substrate, and up to 4 Kgauss applied parallel to the substrate plane. 

Recent Publications: 

1. L. E. Murr, "Structure and Energetics of Vapor-Deposited Indium 
Films," Phys. Stat. Sol. (a),l, 487 (1970). 

2. G. Wong and L. E. Murr, "Structure of Palladium Single -C rys tal 
Films Prepared by Flash Evaporation onto (001) NaCl Substrates, " 
Proc. EMSA, C, J. Arceneaux (ed), Claitor's Publishing Div. 

(Baton Rouge), p. 456 (1970). 

3. K. P. Staudhammer and L. E. Murr, "Thermally Activated Nucleation 
and Faceting in Immiscible Cu-Cr Thin Films, H Ibid, p. 518. 

4. L. E. Murr, "Effects of Substrate Temperature, Pressure, and 
High Evaporation Rates on Nucleation, Epitaxy and Structure of 
Palladium Thin Films, " Thin Solid Films, 7, 101 (1971). 

5. H. P. Singh and L. E. Murr, "Nculeation and Growth Characteristics 
of Palladium and Indium Thin Films, " WESTEC 1971, ASM Technical 
Paper W71-12. 6, 10 pp. (1971). 

6. K. P. Staudhammer and L. E. Murr, "Electron Microscopy of Thin 
Film Immiscible-Miscible Systems Proc . EMSA, C. J. Arceneaux 
(ed. ), Claitor's Publishing Div. (Baton Rouge), p. 218 (1971). 

7. L, E. Murr, O. T. Inal and H. P. Singh, "Nucleation and Epitaxial 
Growth of Vapor-Deposited Films, Electron Transmission and 
Field-Ion Mic roscopy Studies Ibid. , p. 214. 

8. G. I. Wong, H. P. Singh and L. E. Murr, "Electron Diffraction 
Effects in Multi-Layered (Overlapping) Palladium Single-Crystal 
Films , 11 Ibid. , p. 190. 
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9 . 


H. P. Singh and L. E. Murr, "Nucleation and Growth Characteristics 
of Thin Metal Films, ,r Ibid. , p. 216. 


10. L. E. Murr, O. T. Inal and H. P. Singh, "A Field-Ion Microscope 
Study of Vapor-Deposited Platinum, M Thin Solid Films, 9, 241 (1972). 

11. H. P. Singh and L. E. Murr, "Nculeation and Growth Characteristics 
of Palladium and Indium Thin Films, " Met. Trans. , April (1972). 

12. H. P. Singh and L. E. Murr, "Activation Energy and Saturation 
Density for In and Sn Thin Films on NaCl Substrates, " Scripta Met. , 
April (1972). 

13. H. P. Singh and L. E. Murr, "Electron Microscopy of Nucleation 

and G row th of Indium and Tin Films , 11 Phil. Mag. (1972), in the press. 

14. L. E. Murr and H. P. Singh, "Significance of Electric Fields on the 
Growth of Thin Metal Films, " submitted to Appl. Phys. Letters. 

15. H. P. Singh and L. E. Murr, "An Electron Microscope Study of 
Vapor-Liquid-Solid Growth in In and Sn Thin Films, " Proc. EMSA, 

C. J. Arceneaux (ed), Claitor's Publishing Div. (Baton Rouge) (1972). 


1.4.6 Measurement of Interfacial Energies in Solid State 
Metals and Alloys 

N00014-67-A-0269-0010, Office of Naval Research 

L. E. Murr, R. J. Horylev, O. T. Inal, W. N. Lin, 

G. I. Wong 

Objective: 

To clarify points raised concerning the temperature coefficients for 
interfacial free energies in an alloy system by a detailed, systematic 
measurement of the surface, grain boundary, and coherent twin boundary 
free energies in 304 stainless steel over a range of elevated temperatures. 

Recent Progress : 

This work will be completed upon publication of the final report cited 
below. The following abstract of the report describes the work and 
results for the full contract period. 

The method of zero creep of small diameter 304 stainless steel wires 
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coupled with scanning and transmission electron (shadowgraphic) 
microscopy of grain boundary-surface intersection grooves was used to 
measure average values for the surface free energy, Fg, and the grain 
boundary free energy, Ygb* at temperatures of 1060, 1160, 1260 and 
1360°C, from which the temperature coefficients of the specific inter- 
facial free energies, dFg/dT and dygb/^T were measured to be -1.760 
and -0.890 ergs/cm^oC respectively. Mean values of coherent twin 
boundary free energy/grain boundary free energy ratios, Yth^Yab were 
determined for low- torque twin-grain boundary intersections in thin foils 
of 304 stainless steel using techniques of transmis sion and diffraction 
electron microscopy as 0. 01 6, 0.017, 0. 024, and 0. 024 at 910, 985, 

1969 and 1180°C, and values of Ygb substituted at temperatures from the 
measured values of Ygb to obtain average coherent twin boundary free 
energies at temperature. From this data, the temperature coefficient 
of twin boundary free energy was determined to be +0. 007 ergs/cm^°C. 

By assuming the temperature coefficients measured at elevated tempera- 
tures to extend linearly to room temperature, room temperature (2 5°C) 
values of Fg = 3975 ergs/cm^, Ygb ~ 1^70 ergs/cm^, and Y t b = 10 ergs/cm^ 
were obtained; and the room temperature value of stacking-fault free 
energy, Ygp> was determined to be 20 ergs/cm^ by m ultiplying y^b a 
factor of 2. The differences in sign for the temperature coefficients of 
surface and grain boundary free energies compared with the temperature 
coefficient of twin boundary free energy suggests a phenomenologically 
different segregation thermodynamics at these interfaces. 


The techniques of transmission electron microscopy have been 
employed in a study of coherent annealing twins and their energetics in 
ultra-pure aluminum. Coherent twin boundaries have been unambiguously 
identified, and grain-corner twins in aluminum observed for the first 
time. A value of Ytb/ Ygb = °-23 at 4 50°C was obtained for pure aluminum. 

Details of the results are obtained in the final report and publications 
listed therein. In addition, the final report expands and discusses the use 
of techniques developed in this investigation for the measurement of 
interfacial free energies in connection with composite and dispersion- 
hardened systems. Particular emphasis on techniques of scanning 
electron microscopy are presented. 


References and Recent Publications: 

A full list of publications is contained in the following final report 
which can be obtained from the National Technical Information Service or 
the Defense Documentation Center: 


L. E. Murr, "Measurement of Interfacial Energies in Solid State 
Metals and Alloys, " USCEE Report No. 417, July 1972. 
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1.4.7 Influence of Stress on Precipitation and Coarsening 


in Alloys 

AFOSR-71-2087, Air Force Office of Scientific Research 
S. M. Copley 


Objective: 

To achieve a more complete understanding of the effect of stress on solid 
state processes and its concomitant effect on microstructure and proper- 
ties by studying the effect of applied and residual stresses on precipitation 
and coarsening in several alloy systems including: 

(1) titanium-base superconducting alloys 

(2) nickel-base alloys 

(3) Alnico permanent magnet alloys 

(4) unidirectionally solidified eutectic alloys 

Approach: 

To apply a uniaxial stress to alloy specimens as they are cooled through 
their solvus temperature or are aged at a temperature suitable for 
coarsening. Residual stresses are controlled by alloying. The resulting 
microstructures are determined by transmission electron microscopy. 
Particular attention is being given to the possibility of employing "stress 
aging" to produce alloys with novel microstructures and properties. 

Recent Progress: 

In the past year, much progress has been made in achieving the goals 
outlined in our initial proposal. All equipment needed to carry out the 
proposed research has been fabricated or purchased, with the exception 
of the induction melting and crystal growing facility to be completed this 
fall. Arc melted buttons of several Ti-base and nickel-base alloys have 
been prepared. Transmission electron microscopy has been employed 
to determine precipitate morphologies of these alloys aged normally and 
aged under an applied stress. A theoretical study is nearly completed 
that predicts the morphology of coherent precipitates in alloys after 
both normal aging and stress aging, which takes into account such 
factors as lattice misfit, elastic inhomogeneity, elastic anisotropy, 
volume fraction of precipitate, precipitate size and interphase boundary 
energy. A further theoretical study is nearly completed that predicts the 
amount of stress induced segregation that can occur in multi- component 
systems subjected to an applied uniaxial stress. 
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References and Recent Publications: 


1. S. M. Copley, "Application of Continuum Elastic Models to 
Transformations in Solids", presented at the Symposium titled 
"Effect of Stress on Reactions in Solids", 4th Annual Spring Meeting 
of The Metallurgical Society of AIME, Boston, Mass. , May 8, 1972. 

2. G. H. Narayanan and S. M. Copley, "A Composite Model for 
Characterizing Residual Strains in Alloys with Coherent Precipitates", 
in preparation. 

3. M. Menon and S. M. Copley, "Stress Induced Segregation in Binary 
Alloys", in preparation. 

4. S. M. Copley, "The Effects of Interface Coherency on the Morphology 
and Properties of Metallic Composities ”, an invited paper to be 
presented at a Symposium on Interfacial Bonding and Fracture in 
Polymeric, Metallic and Ceramic Composites, November 1972, 
sponsored by the American Chemical Society and U.C.L.A. The 
proceedings will be published in The Journal of Adhesion. A hard 
cover version will be published by Gordon and Breach, Science 
Publishers, Inc. 


1.4.8 Materials Research on High Field Superconductors and 
Rare Earth Cobalt Magnetic Materials 
GH-34652, National Science Foundation 
AT (04-3 )- 1 13, #23, Atomic Energy Commission 

Y. B. Kim, R. Wang, and J. Savage 

Objective: 

To investigate the metallurgical micros tructure in relation to the 
superconducting properties of high field superconductors and the 
coercive forces of rare earth cobalt magnetic materials. 

Approach: 

As for the superconducting materials , we have developed powder sintered 
compacts in which Ng^Sn pow ; der particles are dispersed in copper matrix. 
The ac losses of such superconducting materials may be significantly 
lower than in the conventional multi-filament superconducting cables. 

The metallurgical aspects of these sintered compacts have already been 
described in the last progress report (#14), and currently the measure- 
ments of ac losses in these materials are being carried out. 
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Recent Progress: 


By applying the rapid quenching technique to rare earth cobalt magnetic 
materials, we have obtained coercive forces much larger than attainable 
by the conventional methods. The PrCo 5 and (MM)Co 5 powders obtained 
by splat cooling are found to have very large coercive forces. Small 
grains of 1-5 pi in size are produced in particles of l-40|a. Simple 
binding of these particles under ~27 kOe leads to coercive fields as 
large as 9 kOe for PrCos and 5 kOe for (MM)Co 5 , which are higher than 
the largest values reported for the conventionally fabricated magnets. The 
SmCo 5 powders made by this process, however, do not yield a larger 
coercive field. In splat cooling, particles are subject to much less 
surface defects than that encountered in the conventional grinding 
process. Also, some degree of preferential growth of the grains appears 
to result from the high thermal gradient between the liquid alloy and the 
substrait. In addition, a second phase was found to form surrounding the 
RC 05 phase, which may have improved the surface properties of the RC 05 
phase. Since the presence of the second phase has been shown to have 
advantages to the sintering process for some high quality magnets of 
SmCo 5 and PrCos, it is expected that liquid sintering of fine powders 
prepared by this new process would lead to greater improvement of 
magnetic properties (see Reference a). 

References and Recent Publications: 


1. R. Wang, K. S. Kim, and Y. B. Kim, "Preparation of Fine Powders 
of Rare Earth Cobalt by Rapid Quenching Technique, " presented at 
the 1972 InterMag Conference, Kyoto, Japan. 


1.4.9 Ultr a-Low Temperature Studies of Superconductors 
and Magnetic Alloys 

GP-29031, National Science Foundation 
Y. B. Kim, M. D. Daybell 


Objective: 

I. Understanding of the properties of magnetic alloys from an interacting 
local moment viewpoint. 

II. Investigation of fluctuation dominated properties of type II super- 
conductors. 
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Approach: 

I. Properties of isolated local moments and pairs of local moments in 
a metal are determined from magnetic and transport properties of 
alloys in the temperature range where phonons are frozen out. 

II. Microwave impedance and other transport properties of type II 

superconducting alloys are examined down to very low temperatures 
and up to fields 50 kilogauss. 

Recent Progress: 

Investigation on the interaction between magnetic impurities 
(interaction between Kondo systems) in zero external magnetic fields has 
been described in the previous progress report (#14). Work is continuing 
to extend this investigation in the presence of external magnetic fields. 

In order to resolve the long-standing controversy between the 
Caroli-Maki theory and the Thompson theory, we have measured the 
microwave surface resistance of various Nb-Ta and Pb-In superconducting 
alloys using a derivative technique to determine its rate of change with 
field near Previous dc measurements of this slope have generally 

been interpreted to show agreement with the Caroli-Maki theory, which 
Thompson has found to be incomplete. Measurements of the derivative of 
the normalized surface resistance of samples of dirty type II super- 
conductors with respect to the magnetic field, compare favorably with 
the theoretical predictions of Thompson when the measured values of 
the derivative are extrapolated to H(^2* If the derivative is obtained 
from the region below the transition where its value is relatively constant, 
results which support the Caroli-Maki theory are obtained. Data 
collected without the benefit of the derivative technique would tend to 
yield such results. However, the data obtained by an appropriate 
extrapolation in the immediate vicinity of Hq2 support the Thompson 
correction to the Caroli-Maki theory (Reference a). 

References and Recent Publications : 

1. R. J. Pedersen, Y. B. Kim, and R. S. Thompson, "Microwave 

Surface Resistance of Type II Superconductors in the Vortex State", 
submitted for publication in Physical Review, 1972. 
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1.5 SOLID STATE DEVICES 


1. 5. 1 Soli d State Microwave Devices 

F44620-71 -C -0007, The Joint Services Electronics Program 

H. L. Stover, G. T. Culbertson, B. J. Gordon 


Objective: 

To investigate present devices for new understandings and to investigate 
fabrication and assembly of new devices both leading to better microwave 
devices; to provide service work to other research groups. 

Approach: 

To fabricate microwave devices directly from raw materials using new 
techniques; grow high-quality epitaxial silicon, assemble multiepitaxial 
devices, and study old and new devices over wide frequency ranges. 

Recent Progress : 

Much work was devoted to mastering some of the peripheral processes 
and the frabrication and installation of equipment. All the technologies 
developed in the laboratory over the past several years were merged to 
successfully fabricate and test a Punched-Through Injection Transit-Time 
Microwave Oscillator. Theoretical results have indicated that the device, 
by virtue of its low noise, will make a significant mark on the microwave 
semiconductor device field. Several characteristics of the oscillator are 
such that it could potentially supplant IMPATT and Gunn devices in many 
applications. Several different oscillator structures were fabricated, 
each potentially more efficient. 

The laboratory facilities continue to provide service work for other 
research groups as well as to study new devices and techniques. 

References and Recent Publications: 


1. "A Punch Through Injection Diode: A New Microwave Source and 
Technology Related to Its Fabrication”, G. T. Culbertson, Ph. D. 
Dissertation, 1972; reviewed during the February 1972 JSEP meeting 
at USC. 

2. "I-V Characteristics of Punched-Through Injection Mic rowave 
Oscillator”, B. J. Gordon and H. L. Stover, in preparation. 
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1. 5.2 Microwave Acoustic Devices 


N00014-67-A-0269-0018, Office of Naval Research 
F44620-71 -C-0007, The Joint Services Electronics Program 


K. M. Lakin 


Objective; 

The Microwave Acoustic and Solid State Device Laboratory is directed 
toward a complete investigation of acoustic wave propagation, materials 
growth and evaluation, and a study of acoustoelectric device physics using 
surface elastic waves. 

A pproach; 

By theoretical analysis of wave propagation and interactions with semi- 
conductors and by experiments performed on surface wave delay lines and 
amplifiers . 

Recent Progress; 

Improvements have been made in the computer programs to increase 
the computation speed and accuracy for surface wave boundary value 
problems. An extensive analysis of surface waves on sapphire was carried 
out. The case of ZnO epitaxial films on R plane sapphire was computed 
and the existence of a new high coupling constant layered mode was 
verified quantitatively. This new mode should increase the efficiency of 
signal processing devices and will be examined in greater detail. 

Epitaxial films of AIN and R plane sapphire were grown with improved 
surface finish and quality. This material was used to make a 400 MHz 
surface wave delay line. This high velocity material combination is being 
further evaluated for signal processing devices at 1 GHz where larger 
device bandwidth are obtainable. In addition, silicon- on-sapphire films 
have been grown and research is continuing in an effort to improve the 
semiconducting properties. 
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1. 5. 3 


Ult rahigh Frequency Stable Signal Sources 

F44620- 71 -C -0007, The Joint Services Electronics Program 


K. M. Lakin 


Objective: 

To study a new technique for obtaining stable ultrahigh frequency crystal 
controlled oscillators. The research would involve a study of both 
materials and devices employing surface wave delay lines. 

Approach: 

Surface wave delay lines will be employed as feedback elements in a 
large multiple 2 tt phase shift oscillators. The temperature coefficient 
properties of useful piezoelectrics will be examined both theoretically and 
experimentally to obtain "fundamental mode" oscillations up to 1 GHz. 

Recent Progress : 

This is a new JSEP task proposed on the basis of initial theoretical 
and experimental results at 100 MHz obtained under the broad task 
"Microwave Acoustic Devices" as described in the preceding article. 

Recently we have used a surface wave delay line to phase lock a 
laboratory voltage controlled signal generator as a preliminary step in the 
research. In addition, the phase locked loop has been found useful for 
evaluating the temperature properties of materials for other than 
oscillator applications such as tapped delay lines. 
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1. 6 DEFECT CHEMISTRY AND SOLID STATE ELECTROCHEMISTRY 

1.6.1 High Temperature Defect Chemistry of A l^Q., 

F 33615-72-C-1123, U.S. Air Force (AFSC) 

F. A. Kroger, A. Deganin, B. V. Dutt, and J. Yee 


Objective: 

To determine the defect structure of pure and doped sapphire. 


Approach: 

a) High- temperature conduction studies of AI 2 O 3 , pure and undoped 
with Mg, Co, Ti. 

b) Electromotive force measurements on cells Pt, 0 2 I AI 2 O 3 | Pt , air 
and cells W, H 2 ,H 2 C>| Al^ | W,H 2 ,H 2 0 . 

c) Measurement of absorption and electron spin resonance on A^O^ -Co 
with cobalt in various valency states (2+, 3+). 

Recent Progress: 

Emf measurements on cells Pt, C >2 I AI 2 O 3 I Pt, air have shown that 
the ionic transference number = 1 up to 1450°C. The valency of cobalt 
can be changed from 3+ to 2 + and V. V. by heating in H 2 + HgO and air, 
respectively at 1500°C. Since heating in a low p 02 (e. g. , a CO-CO 2 
mixture) does not have the same result, this indicates that H is entering 
the crystal. 

~ (gas ) + Co 3+ - + (Co^j) 


Presently measurements are being carried and on AI 2 O 3 in H 2 * viz. , 

a) conduction 

b) emf of cells AI 2 O 3 | W,H 2 #H 20 . 

From the latter, both tj = t^ and t^, the proton transference number can 
be determined. 


Present results indicate that both t^l and t^ > 0. At T > 1600°C, the 
valency of Co can be changed in the absence of H 2 . Electron spin resonance 
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O ! 

shows the signal due to Co increasing with decreasing pC> 2 . Optical 
absorption shows a corresponding increase in Co^ + absorption peaks and 
a decrease in Co^ + absorption. Attempts are being made to make the 
measurements quantitative, i, e. , determine oscillator strength. Then 
absorption strength will give [ Co^ + ] and combination with conduction 

<T = C Al 3qU Al with = J c Cq 2 + gives the mobility of 

interstitial Al with a greater accuracy than we have now. 

References and Recent Publications: 


1. R. J. Brook, J. Yee, and F. A. Kroger, "Elec trochemical Cells 

and Electrical Conduction of Pure and Doped AI 2 O 3 , " J. Am. Ceram. 
Soc. 54, 444-451 (1971). 
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2. PLASMAS AND APPLIED ELECTROMAGNETICS 


2.1 Radiation, Waves and Transport Properties in Plasmas 
GK24877, National Science Foundation 


H. H. Kuehl 


Objective: 

To obtain an understanding of the near fields and input impedance of 
radio frequency probes or antennas in warm magnetized plasmas. 

Approach: 

The theory is based on the Boltzmann equation using the electrostatic 
approximation to obtain simple results. Measurements are performed on 
small probes or dipoles in a pulsed r. f. discharge using argon gas. 

Recent Progress: 

A theory has been derived for the input impedance of a short, cylindrical 
dipole antenna immersed in a warm anisotropic plasma which is described 
by the kinetic theory (Boltzmann equation). For a dipole oriented parallel 
to the static ira gnetic field, the input impedance is obtained as a one- 
dimensional integral suitable for numerical integration. Under certain 
conditions , the integral can be evaluated giving analytical formulas valid 
near the second and third electron cyclotron harmonics. The results show 
new contributions to the input impedance due to the excitation of cyclotron 
harmonic waves which propagate near the harmonics of the cyclotron 
frequency. These effects are not predicted by either the cold or hydro- 
dynamic theories. 

Impedance measurements have been performed on a dipole immersed 
in a pulsed, r. f. discharge, argon plasma. The effects of the plasma sheath 
are found to be very important. The experimental results show that the 
input impedance undergoes frequency passes through the harmonics of the 
electron cyclotron Irequency which is in agreement with the predictions of 
the theory. 

References and Recent Publications: 


1. H. H. Kuehl and A. L. Peratt, n T ransmission and Reflection of a Wave 
Obliquely Incident on a Nonuniform Magnetized Plasma", Radio Sci. 7, 
309 (1972). 



2. H. H. Kuehl and A. L. Peratt, "The Plasma Capacitor in a 
Magnetic Field' 1 , Phys. Fluids 15, {1972). 


2 . 2 Integral-Equation Analysis of Axially Symmetric Reflectors 
F44620-71-C -0007, The Joint Services Electronics Program 


W.V.T. Rusch 


Objective: 

Solution of the problem of scattering from axially- symmetric bodies 
using moment techniques to solve the appropriate integral equation for the 
particular problem. 

Approach: 

Application of numerical techniques to invert Maue's Integral Equation for 
the currents on an axially symmetric reflector. 

Recent Progress : 

Maue's integral equation has been formulated for scattering from a 
surface of revolution with a linearly-polarized incident spherical wave. 

The equations have been converted into linear, algebraic matrix equations. 
These matrix equations have been inverted using a computer program 
which has been debugged and checked out using comparisons with classical 
scattering results for a perfectly conducting sphere. The program has 
been applied to scattering from a hyperboloidal subreflector in a 
Cassegrainian feed system. The results do not appear to converge as 
the subreflector thickness vanishes. This is an intrinsic property of an 
H-field integral equation. 

Attention has been focused on the E-field integral equation for the 
same geometry. This equation has a somewhat less diagonally-dominant 
system matrix. However, it is less subject to convergence problems in 
the case of vanishing reflector thickness. This program is also programmed, 
using a program developed by Mautz. The program has been checked out 
for a sphere, but some difficulties are being experienced when applying it 
to a hyperboloid. Work is continuing on this aspect of the problem. 

References and Recent Publications: 


1. "Net-field Polarization in a Magnetical ly Biased Plasma", by 

V.T. Rusch and C. T. Stelzried, has been revised and resubmitted 
to Radioscience. 
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2 . 


"■Applications of Source Multipole Expansions to Simplify Reflector 
Antenna Calculations", by W. V.T. Rusch; presented as a seminar 
at UCLA on May 23, 1972, and to be submitted to URSl/GAP 
Symposium in December, 1972, 

3. "Scattering from a Surface of Revolution Using Integral Equations", 
Engineer's Thesis in preparation by H. A rain. 

4. "Evaluation of the Equivalent Paraboloid Concept", paper in progress 
by W. V.T. Rusch and R. Booth. 

5. "Radiation from a Paraboloid with an Axially Defocused Feed", by 
P. T. Ingerson and W. V.T. Rusch, paper submitted to IEEE 
Transactions on Antennas and Propagation. 
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3. ELECTRICAL ENGINEERING SYSTEMS AND 


INFORMATION SCIENCES 


3. 1 CONTROL SYSTEMS 


3. 1. 1 Basic Studies in Traffic Responsive Control Systems 
GK 24520, National Science Foundation 


N. E. Nahi and H. J. Payne 


Objective: 

To develop methods for the design of freeway traffic surveillance and 
control systems. 

Approach: 

To develop a mathematical model representing traffic flow on freeways 
and use this model to derive optimum control algorithms. Further, 
to derive optimum procedures for estimating traffic variables such as 
average car speeds and car density associated with any section of the road. 

Recent Progress : 

a) Models of Freeway Traffic. We have developed two distinctly 
different types of traffic models which provide complementary tools for our 
research on control and estimation algorithms. The first of these is a 
microscopic model, implemented on a digital computer, in which move- 
ments of individual vehicles are followed. 

The second model developed is an aggregate variable model. In this 
model, the freeway is considered to consist of continuous sections with 
uniform geometric features, each of length approximately 1/2 mile; the 
variables of the model are the density and mean speed of vehicles in each 
section and the flow rate past boundaries of the sections. This model 
provides a capability of simulating an extended segment of freeway for 
several hours with relatively small computational expense. 

b) Regulation of Freeway Traffic . Assuming nominal conditions on a 
freeway, e. g. , those determined from the aggregate variable model using 
fixed ramp metering rates, one can obtain a linearized aggregate variable 
model whose variables are perturbations from nominal values. We have 


- 58 - 



taken as the objective of freeway regulation the maintenance of these 
nominal conditions. This can be achieved by selecting controls (adjustments 
in ramp metering rates) which minimize a cost expressed as a quadratic 
form in the perturbation variables. Suboptimal controls (i.e., those which 
nearly minimize this cost) have been determined with the desirable 
structural feature that each on ramp metering rate is determined by 
traffic conditions in a local segment of the freeway (of length typically 

miles). This design method will be used to develop on ramp control 
software in the experimental Hollywood Freeway project. 

Results of this research have been and will be reported in references. 

c ) Estimation of Traffic Variables. The basic aggregate traffic 
parameters are section mean-speed and section density. These are 
defined as the average of instantaneous speeds of automobiles in a specific 
section of the road, and the number of vehicles per lane-mile in that 
section, respectively. A direct measurement of these variables is 
usually very complicated and clearly not feasible as part of a real time 
traffic control or surveillance procedure. Such measurements, however, 
can be made and may, for instance, take the form of processing successive 
visual images of the road section through aerial photography. The 
complexity of this required processing makes procedures of this type 
infeasible for real time operation. 

The alternative is to perform local measurements at various spatial 
intervals along the road and utilize this information to estimate the 
section mean-speed and section density. The most common local 
measurement device in use today is the so called "pres ence -detector ". 

This device produces an output which is proportional to the time a 
vehicle is located above the location of the sensor (typically a 6‘ x 6' 
square in each lane of a road). 

During the course of the present grant, optimal section mean-speed 
estimators and section density estimators were determined. These 
estimators are designed to utilize the information content of the output of 
the presence detectors located upstream of traffic at the beginning of the 
section. The estimators are of recursive nature and directly operate on 
the output of the presence-detectors. Consequently, they can be 
implemented very simply and in real time. The performance of the 
estimators was evaluated using (1) a car following simulation and 
(2) actual data derived from aerial photographs. The results of this 
study appear in detail in reference, a copy of which is attached to this 
proposal. The results, both from the point of view of simplicity of 
the estimators and quality of the estimates are very encouraging. 

References and Recent Publications: 


1. Isaksen, L. , and H. J. Payne, "Regulation of Freeway Traffic", 
submitted to IEEE Trans. Automotive Control. 
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2. Payne, H. J. , "Models of Freeway Traffic and Control", Simulation 
Council Proceedings, Vol. 1 (Mathematical Models of Public Systems), 
Chapter 6. 

3. Mikhalkin, B. , H. J. Payne, and L. Isaksen, "Estimation of Speed 
from Presence Detectors", to appear in the Highway Research Board 
Record, Freeway Operations. 

4. Nahi, N. E. , and Trivedi, A., "Recursive Estimation of Section 
Mean Speed and Density", USCEE Report No. 414, The Electronics 
Sciences Laboratory, University of Southern California (1972). 


3.1.2 Studies of Large Scale Linear Systems 

F44620-71 -C -0007, The Joint Services Electronics Program 


H. J. Payne 


Objective: 

To develop techniques to take advantage of invariance properties of 
large-scale linear systems. 

Approach: 

The studies will relate invariance properties of the underlying physical 
system to exploitable algebraic structure of the linear system model. 

Recent Progress: 

Linear dynamic systems with a periodic structure can be described 
by equations of the form 


x = Ax + Bu 
y = Cx 

in which A, BB' and C’C are circulant matrices. The optimal regulation 
of such systems is then specified by 

u = B 1 Px 

in which P is a circulant solution of the algebraic Riccati equation 

PA + A'P - PBB'P + C'C = 0 

These notions have been generalized to include matrices consisting of 
blocks of circulants. 
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The circulant properties can be exploited to yield a set of low order 
Riccati equations whose solutions provide the desired optimal regulator. 
This reduction in order has been shown to result in considerable 
computational savings in solving the algebraic Riccati equation. 

This work was reported in the reference cited below and is considered 
to be complete. Extensions of the methodology employed here to 
generalized Toeplitz systems are being considered. 

References and Recent Publications : 

1. H. J. Payne, "Circulant Solutions of the Riccati Equation", submitted 
to IEEE Control Systems Society Transactions. 


3.1.3 Input- Output Structure of Linear Systems 

DA-ARO-D-3 1- 124-71-G48, U.S. Army Research 
Office - Durham 


L. M. Silverman and B. C. Moore 


Objective: 

The main objective of this project is to characterize the input-output 
mapping properties of multi-input, multi-output linear systems and to 
apply these characterizations to various synthesis problems of control, 
estimation coding, etc., which involve systems which can be modeled 
appropriately. 

Approach: 

The basic approach is to utilize the structure algorithm developed by the 
principal investigator as a means of characterizing various quantities 
associated with a linear system. This has proved successful in attacking 
such problems as system inversion, input-output decoupling and model 
ma tching. 

Recent Progress : 

During the reporting period work was completed on the general model 
matching problem [l]. A necessary and sufficient condition was found 
for determining when an existing system can be dynamically compensated 
so that it has the same input-output response as some desired model and a 
procedure based on the structure algorithm for designing the required 
compensation was given. One potential application area for these results 
is in flight simulation where it is desired to compensate an existing 
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aircraft so that it flies like an as yet unbuilt one. 

Another problem presently being considered is that of characterizing 
the ‘'zeros” of a multi-input multi-output system. Such "characterizing 
atoms are important in many system problems. For the class of systems 
which describe convolutional codes a complete solution has been given [2], 
Work is nearing completion for the general case. 

References and Recent Publications: 

1. B. C. Moore and L. M. Silverman, "Exact Model Matching by State 
Feedback and Dynamic Compensation", presented at the 1971 Decision 
and Control Conference, December 1971. Also to appear in the IEEE 
Transactions on Automatic Control. 

2. B. C. Moore and L. M. Silverman, "A New Characterization of 
Feedforward, Delay-Free Inverses ", pres ented at the 1972 Inter- 
national Symposium on Information Theory, Pacific Grove, Calif. , 
January 1972, 


3.1.4 Optimal Control Problems with Mixed Control - Phase 
Variable Equality and Inequality Constraints 
A USC Project - No Contract Support 

L. W. Neustadt 

Objec tive; 

To obtain necessary conditions of maximum- principle type for optimal 
control problems containing equality and inequality constraints involving 
both the control and phase variables. 

Approach: 

The problem was attacked from the viewpoint of a mathematical program- 
ming problem in Banach spaces with an equality constraint in an infinite 
dimensional Banach Space (La,). 

Recent Progress: 

Necessary conditions in the form of a maximum principle of the 
Pontryagin type were obtained for an optimal control problem containing 
ordinary differential equation constraints of the usual form 

x(t) = f(x(t), u(t), t), t 1 < t < t 2 , 


- 62 - 



phase inequality constraints of the form 


p l (x(t), u(t), t) <0 , i = tel 1 , 

p l (x(t), u(t), t) = 0 , i = (j’+l, . . . , |j , tel', 

where I 1 is a given subset of [ t, t^ J. In addition, more conventional 
constraints of inequality type or the control variable 

q 1 {u(t), t) <0, i = 1, . . . , o , 

as well as constraints on the initial and final values of x were also included. 

The results obtained were similar to previous ones of other authors, 
but much weaker "compatibility” conditions on the qi and the pi were 
assumed in this study than was the case in the earlier work. 

References and Recent Publications : 

1. 1 paper completed and submitted to SIAM Journal on Control: 

"Optimal Control Problems with Mixed Control- Phas e Variable 
Equality and Inequality Constraints", by K. Makowski and 

L. W. Neustadt. 

2. Book - completed and submitted to publisher (Holt, Rinehart, & 
Winston): Optimization: A theory of Necessary Conditions, by 
L. W. Neustadt. 


3.1.5 Simplified Filter Models 

F44620-71-C-0007, The Joint Services Electronics Prograim 


D. D. Sworder 


Objective: 

To derive algorithms which will provide near optimal filters with reduced 
computational requirements. 

Approach: 

The filters will be derived using projection methods developed in the 
earlier reporting period. 
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Recent Progress : 

An algorithm has been derived for finding the best estimate of a linear 
function of the state variables of a linear system when a constraint is 
placed on the computational complexity of the estimator. Necessary 
conditions on the estimator dynamics have been derived. When the dynamic 
modes of the estimator are prespecified, a constructive method for finding 
the bewt filter has been given. Without this restriction a gradient search 
must be made in conjunction with the aforementioned method to find the 
true optimum. 

References and Recent Publications: 

1. "Optimal Sub-optimal Linear Filtering 1 ', C. B. Shortwell, Ph. D. 
Dissertation, 1972. 

2. "Sub optimal Control of Linear Systems Derived from Models of Lower 
Order", 1 97 1 Joint Automatic Control Conference, pp. 926-935, 1971; 
also AI A A Journal, Vol. 9, No. 8, pp. 1461-1467, August 1971 (with 
R. O. Rogers). 

3. "Constrained Linear Filtering as a Multistage Process", 1971 Joint 
Automatic Control Conference, pp. 135-139, 1971 (with A. Jain). 

4. "Feedback Estimation Systems and the Separation Principle of 
Stochastic Control", 1971 Joint Automatic Control Conference, pp. 834 
840, 1971; also IEEE Trans . on Automatic Control, Vol. AC-16, 

No. 4, pp. 350-354, August 1971 (with S. C. Sworder). 

5. "Optimal Control, Repair and Inventory Strategies for a Linear 
Stochastic System”, 1971 IEEE Conf. on Decision and Control, 
pp. 393-394, 1971 (abstract); also IEEE T rans . on Systems, Man 
and Cybernetics, July 1972 (with T. Kazangey). 
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3.2 COMMUNICATION SYSTEMS 


3.2.1 Nonparametric and Adaptive Techniques in Communication 
GK 14190, National Science Foundation 

NAS 5-21505, National Aeronautics and Space Administration 

NGL-05-0 18- 1 18, National Aeronautics and Space 
Adminis tration 


L. D. Davisson 

Objective: 

To investigate and evaluate the performance of nonparametric and 
adaptive techniques in communication applications; to develop a practical 
data compression system for video data. 

Recent Progress : 

The effect of small sample sizes on rank tests has been investigated. 

It is found that rank tests perform nearly as well as the optimum detector 
for charge signals in Gaussian noise. 

A continuous time RC filter-limite r-RC filter system has been 
examined. An important conditional mean linearity conjecture has been 
disproven. 

The effect of dependence on adaptive filter performance has been 
partially evaluated when the dependence is confined to a finite time interval. 

A real time data compression system has been developed and put 
into partial operational use for the transmission of video data between 
Wallops Island, Virginia, and Suitland, Maryland. 

References and Recent Publications : 

1. L. D. Davisson and R. L. Kutz, "A Real Time Programmable Data 
Compression System", Proc. Int'l Telemetering Conference, 1971. 

2. L. D. Davisson, "A Survey of Data Compression", Proc. Allerton 
Conf. on Circuit and System Theo. , 1971. 
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3.2.2 Graph Theory, Networks, and Mathematical Modeling 
with Systems Engineering Applica ti ons 

AFOSR 68-1555, The Air Force Office of Scientific Research 


S. W. Golomb, I. S. Reed, H. Laiten and R. A. Scholtz 
Objective and Approach: 

To investigate the many new and continuing problems which are out- 
growths of previous work under the general title of "Automata and 
Communications" with specific attention directed to three topics listed 
in the following section. 

Recent Progress: 

a. Coding Problems in Information Theory. Publications and 
lectures on this topic for the past three years include: 

1) "Perfect Codes in the Lee Metric' and the Packing of Polyominoes, " 

by S. W. Golomb and L. R. Welch, SIAM Journal on Applied Math, 
vol. 18, no. 2, January, 1970, pp. 302-317. 

2) "Algebraic Coding and the Lee Metric, " by S. W. Golomb and 
L. R. Welch, Chapter 10 of Error Correcting Codes , edited by 
Henry B. Mann, John Wiley and Sons , 1969, pp. 175-194. 

3) "A General Formulation of Error Metrics", by S. W. Golomb, IEEE 
Transactions on Information Theory, vol. IT-15, no. 3, May 1969, 
pp. 425-426. 

4) "Tiling with Sets of Polyominoes," by S. W. Golomb, The Journal 
of Combinatorial Theory, vol. 9, no. 1, July 1970. 

5) "Sphere Packing, Coding Metrics, and Chess Puzzles," by 

S. W. Golomb, presented to the Second Chapel Hill Conference on 
Combinatorial Mathematics and Its Applications, at the University 
of North Carolina, May 1970, and published in the Symposium 
Proceedings, August 1970. 

6) "Combinatorial and Geometric Aspects of Coding, " Research 
Award Address by S. W. Golomb to the University of Southern 
California Faculty, May 8, 1969. 

7) "Geometry of Coding," invited address by S. W. Golomb to the 
National Telemetry Conference, April 27, 1970, Los Angeles, 
California. 
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8) ’’Information Theory and Coding, " invited hour-address by S. W. 
Golomb to the SIAM meeting at Caltech, 24 April 1971. 

9) ’’Efficient Coding and the Packing of Octahedra, " invited talk by 
S. W. Golomb at the Seminar on Structure, USC, 4 May 1971. 

10) "The Systematic Selection of Cyclically Equivalent Codes, " by 

I. S. Reed and C. T. Wolverton, IEEE Transactions on Information 
Theory, March 1972, 

11) "Notes on the Arithmetic BN Modulo A Codes, " by I. S, Reed and 
A.CoL. Chiang, IEEE Transactions on Information Theory, to 
appear, 1972. 

b. Analytical Foundations of Statistical Communication Theory. The 
articles contributing to the mathematical foundations of information theory 
and communication theory starting in 1970 are: 

1) "The Lambda Method in Prime Number Theory," Journal of Number 

Theory, V. 2, N. 2, May 1970, pp. 193-198. 

2) "A Class of Probability Distributions on the Integers, " Journal of 
Number Theory, V. 2, N. 2, May 1970, pp. 189-192. 

3) "Powerful Numbers, " American Mathematical Monthly, vol. 77, 
no. 8, October 1970, pp. 848-852. 

4) "A New Arithmetic Function of Combinatorial Significance, " 

S. W. Golomb, accepted by the Journal of Number Theory, 1971. 

5) "On the Survival of Sequence Information in Filters, " by S. W. 
Golomb, IEEE Transactions on Information Theory, March 1972. 

6) "An Information-Theoretic Measure of Image Similarity, " by 

I. S. Reed and L. A. Liporace, to appear, IEEE Transactions on 
Information Theory, 1972. 

c. Automata, Logical Design and Graph Theory . Publications on 
this topic since 1970 are: 

1) "On the Enumeration of Finite State Synchronous Sequential 

Machines," USCEE Report 389, by I. P. Bottlik (based on Ph. D. 
dissertation supervised by I. S. Reed), June 1970. 
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2) "Error Tolerant Sequential Circuit Design Using Error Correcting 
Codes," USCEE Report 388, by R. W. Larsen (based on Ph.D. 
dissertation supervised by I. S. Reed), June, 1970. 

3) "The Stability of Counting Sequences Under Stage Delays, " by 

S. W. Golomb, jointly with M. Cohn, S. Even, and A . Lempel, The 
SIAM Journal on Applied Mathematics, vol. 20, no. 2, March 1971, 
pp. 183-188. 

4) "K-th Order Near Orthogonal Codes, " by I. S. Reed, IEEE Trans- 
actions on Information Theory, vol. IT-17, No. 1. 

5) "Some Decompositions of the Integers from 0 to p n ~^, 11 by 

S. W. Golomb, American Mathematical Monthly, February, 1972. 

6) "Numbering the Nodes of a Graph, " by S. W. Golomb, an invited 
chapter -length contribution to the book. Computing and Graph Theory , 
Academic Press, edited by R. C. Read, 1972. 

7) "Graceful Graphs, " Scientific American, March 1972, by S. W. 
Golomb. 

8) "The Inner Product Computer, " by E. E. Schwarzlander , Jr., USC 
Ph.D. dissertation supervised by I. S. Reed, June 1972. 

9) "Coding Techniques for Failure -Tolerant Counters, " by I. S. Reed 
and A.C.L. Chiang, IEEE Transactions on Computers, November 
1970. 

10) "Path Sensitization, Partial Boolean Difference, and Automated 
Fault Diagnosis, " by I. S. Reed, A.C.L. Chiang, and A . V. Barnes, 
IEEE Trans, on Computers, February 1972. 

11) "Redundancy by Coding Versus Redundancy for Failure-Tolerant 
Circuits," by I. S. Reed and R. W. Larson, IEEE Transactions on 
Computers, February 1972. 

12) "Symbolic Design Techniques Applied to a Stored Microprogram 
Computer , " by I. S. Reed, COMPUTER, May/ June, 1972. 

13) "Boolean Difference Calculus and Fault-Finding, " by I. S. Reed, 
to appear, SIAM Journal. 

14) "On-Off Digital Communications," by I. S. Reed and R. A. Scholtz, 
presented at Asilomar Conference (Jaunary 1972) and submitted 
for publication. 
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3.2.3 Frequenc y Tracking Study 

F44620-71 -C -0007, The Joint Services Electronics Program 
DA -ARO-D-3 1-124, Army Research Office, Durham 


C. L. Weber and R. L. Scholtz 


Objective: 

To determine the performance of approximate optimal estimators of 
frequency from the observation of a frequency modulated, fading signal 
in wideband noise. 

Approach: 

Use nonlinear recursive estimation theory to determine the estimator 
and simulator to check the real performance against performance of a 
linearized model. 

Recent Progress : 

This project is now completed with results published as reported 
below. 

However, as an off-shoot of this study, examination was made of 
the performance of an on-off communications system, which required 
only frequency synchronization, i. e. , not phase synchronization, in 
its mechanization. The results were reported at the IEEE International 
Symposium on Information Theory, Pacific Grove, Calif. (Jan. 1972) 
and submitted to IEEE Transactions on Communication Technology. 

References and Recent Publications : 

1. Dissertation: Thelma Rodriguez, "Frequency Tracking Loop 

Study", Ph. D. , 1971. 

2. T. M. Rodriguez, R. A. Scholtz, C. L. Weber, "Frequency 
Tracking Loop Study", Hawaii International Conference on Systems 
Sciences, Jan. 1972. 
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3.2.4 Signal Acquisition, and Synchronization Studies 

F44620- 7 l-C-0007 , The Joint Services Electronics Program 
N000 14, 67A-0269, Office of Naval Research 


W, C. Lindsey 


Objective : 

To further develop the theory of signal acquisition and synchronization in 
the presence of noise and interference. 

Approach: 

To develop a signal acquisition theory including acquisition time, acquisi- 
tion range and probability of acquisition in the time interval [t Q , t] as a 
function of the system phase detector characteristic, loop filter and signal 
wave forms. 

Recent Progress: 

The transition probability density function, which statistically characterizes 
the phase error process in a first-order synchronization system, has been 
developed using the theory of continuous Markov processes. This statisti- 
cal characterization has applications to establishing the fundamental 
limitations on performance of a sync system operating in noise. More 
recently, work has proceeded on the development of a theory which 
specifies the time to acquire a signal when the synchronization system 
has an arbitrary phase detector characteristic. Emphasis is placed upon 
second and third order sync systems. Some experimental work is 
expected when the new communication laboratory is established. 

Applications of the theory are applicable to design and planning of advanced 
communication, identification, navigation and command. 

References and Recent Publications : 

1. Book entitled. Synchronization Sys terns , W. C. Lindsey, manuscript 
completed, to be published June 1972. 
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3.2.5 Digital Image Enhancement Techniques 
GK-2058, National Science Foundation 


H. C. Andrews 


Objective: 

This study of digital transform image enhancement techniques has 3 

main goals: 

a. An investigation as to the applicability of traditional 
analytic techniques resulting in orthogonal transformations 
to digital image enhancement. This would include super- 
resolution and Karhunen-Loeve principles for the deter- 
ministic and stochastic cases respectively. 

b. An investigation into the possibility of efficient implemen- 
tation of general orthogonal transforms as suboptimum 
solutions to number 1 above. 

c. An investigation into the use of known, efficiently 
implementable transforms for image enhancement. Such 
transforms include the Fourier, Hadamard, generalized 
Walsh, Haar, generalized Haar, and a variety of other 
unnamed transforms, all implementable on the order of 

pN logp N operations or less compared to the N 2 operations 
normally required. Here p n = N. 

Approach: 

a. Catalog characteristics and properties of existing transforms 
suitable for image enhancement. 

b. Investigate the structure of transformed images for specific classes 
of degraded images, with specific emphasis on low contrast images, 
using various transforms. 

c. Evaluate the image enhancement attainable by linear and nonlinear 
processing of transform samples for various transforms. 

d. Investigate potential algorithms and computational techniques for 
those image enhancement transforms with most significant promise. 

Recent Progress; 

Activity during the reporting period was directed to 2 topics: 

a * Image Processing Using Digital Computers. Recently digital 

computers have been brought to bear on the problem of processing 

natural quality imagery with surprising success. Optical processes, 
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communication systems, and automatic image detection are three areas 
of application for the computer techniques. Yet, perhaps the most 
exciting areas of research are those in which image restoration, 
enhancement and measurement are implemented through the power of 
the large scale general purpose digital computer. Details of our 
research can be found in the paper (7, below) designed to present a few 
of these areas of application in the hopes that the reader can get a general 
feel for the flavor and possibilities of combining the computer and natural 
images for processing purposes. Examples might include bandwidth 
reduction, detection, biomedical processing and pseudocolor. Probab 
the technique of digitally generated pseudocolor is in its most formative 
stages and applications are developed to illustrate a variety of examples 
from this relatively new field. 

b. Nonlinear Intrin sic Dimensionality Computations. In the feature 
selection aspect of mathematical pattern recognition, the raw data and 
dimensionality of the measure space is usually very large. Therefore, 
some form of dimensionality reduction has been commonly considered 
as a practical pre-processing method. Based on one of previous methods, 
which increases the variance while maintaining the local structure a 
technique is developed to determine intrinsic dimensional structure. 

The minimal spanning tree in the cluster is studied in order to tie it 
to the cluster structure. It has been proven that if the rank order is 
always maintained, a minimal spanning tree will not change its structure 
in the unfolding process. A cost function is introduced to guide the 
maintenance of the rank order. Two criteria of using the cost function 
to increase the variance have been introduced. The local structure is 
related to the connection on the minimal spanning tree. Several methods 
to define the local region are suggested. A program is implemented and 
tested to find the intrinsic dimensionality of the clusters. To date the 
technique has successfully unfolded helices , hemispheres, cubes, and 
lower dimensional surfaces imbedded in higher dimensional spaces. 

References and Recent Publications : 

1. "Data Compression and Enhancement of Sampled Images", with 
A. G. Tescher, A pplied Optics , Vol. 11, No. 4, pp. 919-925 
(April 1972). 

2. "An Experimental Comparison for Karhunen -Loeve and Other 
Unitary Transforms for Feature Selection Using RDT 1 ', Hawaii 
International Conference on Systems Sciences (January 1972). 

3. "Two-Dimensional Sequency Filters ", with H. F. Harmuth and 
K. Shibata, IEEE Transactions on Communication Technology, 

June (1972). 
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4 . 


"Image Display Enhancement and Maximum Entropy 11 , with 
A. G. Tescher, Hawaii International Conference on System Science, 
January (1972). 

5. "Performance Measures for Transform Data Coding", with W. K. 

Pratt and J. Pearl, IEEE Transactions on Communication Tech- 
nology, June 1972. 

6. "Digital Computers and Image Processing", Invited Paper, Endeavour, 
May (1972). 

7. "Image Processing Using Digital Computers", (with H. G. Tescher, 
and R. P. Kruger), IEEE Spectrum (July 1972). 
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3.3 COMPUTER SYSTEMS AND AUTOMATA THEORY 


3.3.1 Automata and Formal Language Theory 
GJ-28787, National Science Foundation 

S. Ginsburg, A. Gabrielian, D. Kiel, J. Outwater 
and B. Rovan 


Objective; 

(1) To develop suitable models of computer devices and software, and 

(2) to use the mathematical models to answer questions of interest to 
computer science. 

Approach: 

The method consists of studying various families of devices and formal 
languages, both in the specific and in the abstract. 

Recent Progress : 

In the work leading to the report of Reference 1, acceptors with 
nonlinear types of storage were studied. In particular, a model was intro- 
duced which is braod enough to encompass multitape AFA as well as devices 
with more exotic, but practical storage structure, such as trees, linked 
lists, two dimensional storage, etc. By appropriate recoding, it is shown 
that such acceptors still define AFL. 

In preparing the report of reference 2, a class of grammars, called 
context- limited, which are weakly equivalent to the contest free, were 
studied. These grammars may be viewed as the ones for which the scan 
limited automata are the parsing algorithm. It Is shown that the 
unambiguous context limited grammar languages and the unambiguous 
context-free languages coincide. 

During the report period, five papers (3-7) were published. 

References and Recent Publications: 


1. A. Gabrielian and S. Ginsburg, "Structural Storage AFA 1 ', submitted 
for publication. 

2. T. N. Hibbard, "Context Limited Grammars", 

3. S. Ginsburg and J. Hopcroft, "Images of AFL Under Certain 
Families of Homomorphisms", Math, Systems Theory, Vol. 5 (1971), 
pp. 216-227. 
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4. R. Book and S. Ginsburg, ’’Multi Stack Counter Languages”, Math. 
Systems Theory, Vol. 6 (1972), pp. 37-68. 

5. S. Greibach and S. Ginsburg, ’’Multitape AFA ", Journal of A ssoc , 
for Comp. Mach., Vol. 19 (1972), pp. 193-221. 

6. J. Golds tine, "Substitution and Bounded Language”, Journal of 
Computer and System Sciences, Vol. 6 (1972), pp. 9-29. 

7. J. Guiliano, "Writing Stuck Acceptors”, Journal of Computer and 
System Sciences, V 6 (1971), pp. 168-204. 


3.3.2 Decomposition of Scenes Into Bodies Using Distance 
Functions 

F44620-71-C-0007, The Joint Services Electronics Program 


A. L. Zobrist and W. L. Thompson 


Objective: 

To develop algorithms or computer programs which calculate a "distance 
function” between parts of an image or picture which approximates human 
perception of distance. Gestalt clustering of parts into wholes can be 
performed using this distance function. This capability is needed for 
image processing and robot vision applications. 

A pproach: 

Interactive computer graphics will allow rapid development and testing of 
programmed distance functions. Individual photoprocessing will allow 
accurate comparison with human perception. Application to real world 
scenes will follow. 

Recent Progress : 

In this new project, recent activities included the following. The 
interactive graphics program has been written and debugged. Seven 
elementary distance functions have been programmed. Informal psycho- 
physical tests have been performed using linear combinations of these 
distance functions, and results show that this approach will work on 
simple scenes. One interesting result is that the orientation of small 
edges is a strong textural cue for the decomposition of scenes into bodies. 
A graduate student is presently adapting this work for application to real 
world scenes. 
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3.3.3 Hardware Implementation of Operating Systems 

F44620-71-C-0007, The Joint Services Electronics Program 


R. Card 


Objective: 

To study the implementation in hardware of the functions of the operating 
system of a computer. 

Approach: 

Microprogramming will be used to achieve the desired hardware control. 
Recent Progress : 

Various tasks of an operating system have been studied with the 
object of determining how they might best be implemented or augmented 
by microprogrammed routines and what the advantages are of micro- 
programming over conventional implementation techniques. These tasks 
include scheduling and allocation of resources, file operations, inter- 
communication and control of application programs, and interrupt service 
routines . 

Microprogrammed routines have been designed to implement the 
control and intercommunication of parallel processes. The design is 
based on theoretical studies that have appeared in recent publications by 
several authors. The usefulness of the microprogram implementation of 
an algorithm for parallel process control and communication is in the 
reduction in the amount of main memory used by the operating system 
and the increase in the speed of execution of the algorithm. 

In order to further study the performance of microprogrammed 
routines designed in this research, a microprogram simulator has been 
installed on the university's 3 70/155 computer system. Further 
research will continue under other contracts with the goal of studying, 
under simulated execution, the operation of designs for the micro- 
programmed implementation of various tasks of an operating system. 
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3.3.4 Computer Program Translation 

N00014- 67A-0269-0020, Office of Naval Research 


P. Gilbert 


Obj ective: 

To obtain formal models of syntactic analysis and translation which 
include parsing schemes as explicit components which are adequate {in 
a strict formal sense) for the parsing and translation of actual programming 
languages, and which can be "carried into practice" to yield rigorously 
based bur efficient translators. 

Approach: 

To extend previous work wherein models called the analytic grammar and 
the analytic transducer were developed and their properties were 
studied. 

Recent Progress : 

Technical Report No. 1 obtains classes of grammars which are 
adequate for the parsing, and classes of transducers which are adequate 
for the translation of actual programming languages. Subsequent work 
has emphasized the "carrying into practice" of these models. 

Specifically, a construction has been developed for a Turing machine 
parser. This construction exactly parallels, in certain important 
respects, the construction of syntax tables (or segments of program) of 
actual parsers. This construction has been shown to be especially 
amenable to minimization. Measures of efficiency have been outlined 
for this parser, and related to the efficiencies of actual parsers. Some 
simple minimization techniques have been outlined. 

References and Recent Publications : 

1. Technical Report No. 1: "Syntactic Analysis and Translation of 

Programming Languages, " December 1971. (Revised version also 
submitted for publication to the Journal of the ACM) 
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3.3.5 Sy nthesis of Sequential Machines 

F44620-71 -C-0007, The Joint Services Electronics Program 

R. S. Kashef 


Objec tive: 

To provide a set of functions which can synthesize a large class of 
sequential machines that are of minimal complexity and can be easily 
fault tested and corrected. 


A pproach: 


To consider the synthesis problem as a class of four sets of transforma- 
tions. T , T , T and T , with each element of the class consisting 
(*) (2) (3) (4) 


of t... X t._. X t... Xt... where t . is an element of T... 

(1) (2) (3) (4) (t) (i. 


and T ^ transforms 


a given problem into the equivalent minimum flow table, T provides the 
state equations for the given problem and T provides the' minimum cost 
equations based on the given technology. 


Recent Progress : 

Significant results have been accomplished toward the evaluation of 
T( 2 )» with the immediate result being that APC codes can be used to code 
any flow tables up to 1 6- rows of sequential machine and n-rows of flow- 
tables with maximal partions. These codes are n + log 2 n. Toward the 
overall synthesis problem, synthesis in cellular spaces are being studied 
with partial results indicating that Bank's cellular space seems to be 
universal and Shoup's cellular space is not practical for synthesis. New 
spaces are being investigated both for universality and practical simplicity, 
specifically those made up of functions of 2 variables. Also, synthesis 
of sequential machines in cellular space of less complexity (fewer states) 
are being investigated. 

References and Recent Publications: 


1. "A Formal Approach to the State Assignment of Asynchronous 
Sequential Machines”, 1EEETC, Computer Group Repository 
R -71-159, Nov. 1971. 

2. "On the Race-Free Assignments and Classifications of 9 to 16 Rows 
Asynchronous Sequential Machines”, IEEE TC Computer Group 
Repos itory R -71 - 1 58, Nov. 1971. 

3. "On the Race-Free Assignment and Classification of Asynchronous 
Sequential Machines”, Proceedings of the 5th Hawaii International 
Conference on Information and System Sciences, Jan. 1971. 
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4 . 


"Augmented Parity Check Codes for Encoding of Asynchronous 
Sequential Machines", Revised Nov. 1971 for IEEETEC, accepted 
for publication. 


3.3.6 Generation of Fault Detection Tests for Intermittent 
Faults in Sequential Circuits 

F44620-7 l-C-0067, The Joint Se rvices Electronics Program 
N00014-67-A, 0269-001 9, Office of Naval Research 
GK-23886, National Sci ence Foundation 

M. A. Breuer 

Objective: 

a) To determine models for intermittent faults in digital circuits; 

b) To determine an effective procedure for generating tests for 
detecting these faults; 

c) To analyze the statistical properties of detection of intermittent 
faults. 

Approach: 

To model the time domain into a space dimension. Single intermittent 
faults in the time domain (faults that come and go) are thus mapped into 
multiple solid faults. Models will then be developed for pattern sensi- 
tive faults which occur in devices such as dynamic RAM memories. 

Models for delay induced faults will also be developed. 

Recent Progress : 

The appearance of an intermittent failure is often due to reasons other 
than signal line values. For these failures a statistical model independ- 
ent of circuit input and state can be used. In this research, one such 
model was introduced based upon Markov chains. Some of the statistical 
properties of detecting intermittent faults were then investigated. A 
concept of a fault pattern was introduced for characterizing when a fault 
pattern was introduced for characterizing when a fault is present, and 
in which permanent faults and the fault free situation are included in 
special cases. Hence, the model is a significant generalization of 
previous work. We have shown that, in general, no single test sequence 
will detect all possible fault patterns. We have shown how to calculate 
the probability of detection of a test sequence. We indicated a procedure, 
based upon the d-algorithm, for constructing a test for an intermittent 


- 79 - 



failure. We show that the behavioral effect of an intermittent failure 

on a circuit may be to increase the number of states in the circuit. 

References and Recent Publications : 

1. ’’Generation of Fault Detection Tests for Sequential Circuits' 1 , 

M. A. Breuer, Digest of 1971 International Symposium on Fault- 
Tolerant Computing, March 1-3, 1971, pp. 18-21. 

2. "An Algorithm for Generating a Fault Detection Test for a Class 
of Sequential Circuits", M. A. Breuer, Theory of Machines and 
Computations, Academic Press, 1971. 

3. "A Random and an Algorithmic Technique for Fault Detection Test 
Generation for Sequential Circuits", M. A. Breuer, IEEE Trans, 
on Computers, V C-20, Nil, pp. 1364-1370, November 1971. 



3.4 ELECTRIC POWER SYSTEMS 


3.4.1 Intermediate Term Load Forecasting 

P. O. 20238, Los Angeles Department of Water and Power 


J. Lamont 


Objective: 

Investigate presently accepted techniques for load forecasting and report 
which would have been most accurate for the Dept, of Water and Power. 
Determine the significance of various factors and make a final forecast. 

Approach: 

Use extrapolation and regression techniques to forecast future electric 
loads. 

Recent Progress : 

Data has been gathered relating to economic, geographic and climatic 
conditions in Los Angeles. Studies are now in process to determine how 
these conditions effect the electric demand. 
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4. MATHEMATICAL BIOSCIENCES AND BIOMEDICAL ENGINEERING 


4. 1 Mathematics in Biology and Medicine 

GM 16197-05, Department of Health, Education, and Welfare 
GP 29049, National Science Foundation 
AT (04-3)-113, Atomic Energy Commission 

R. Bellman, E. Angel, S. P. Azen, J. Buell, A. Jain, 

R. Kalaba, B. G. Kashef, A. Jain, C. Kell, M. J. Merritt, 

S. Osaki, S. Ueno, R. Vasudevan 

Objective: 

In engineering, mathematics has repeatedly furnished the conceptual 
and analytic mechanisms for the design and effective operation of major 
systems. It has been clear for some time that we must pursue a similar 
program in the study of biological and medical phenomena. The human 
body has the ability to maintain homeostasis, to perform complex actions 
like walking and talking, seeing and hearing, to operate a brain which 
stores and retrieves data and makes decisions, and finally to reproduce 
itself. All of these activities represent feats of astounding dimension. 
Most defy complete explanation at the present time despite the enormous 
effort devoted to their study. Even partial understanding, however, will 
have wide-ranging consequences with rewards to society which will dwarf 
all previous contributions of science. 

A basic part of the challenging set of problems that faces us then is 
that of combining modern technology, as represented for example, by 
computers and radiation physics, with modern mathematics to study the 
multiple questions involved in the prevention and cure of illness, in 
patient care, and in the education and training of medical and paramedical 
personnel as well as research personnel. 

Approach: 

In the many research papers produced in this program, we have 
repeatedly indicated how even simplified versions of medical processes 
lead quickly to quite intricate mathematical equations. Conversely, we 
have described how the abstruse mathematical problems we have treated 
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arise from specific diagnostic and therapeutic questions, and have, in 
many cases, immediate application. 

Initially, it is perhaps surprising to people that the human body 
should produce mathematical complications which so often greatly surpass 
in difficulty the familiar glamorous problems of the world of physics and 
engineering. Because of familiarity with body processes we tend to 
ignore their complexity, and assume that physiological processes, if not 
actually simple, are at least capable of simple treatment using familiar 
methods. Unfortunately, this is not the case. 

It is only when we attempt to produce the numbers required for the 
many decisions in medical treatment that we begin to appreciate the true 
complexity of the processes of biology and medicine. 

Recent Progress: 

The following paragraphs provide a survey of the broad areas of 
interest to this research group and the highlights of recent accomplishments. 

a. Pharmacokinetics and Chemotherap y 

As the result of pioneering research by a number of distinguished 
physiologists, Teorell, Kruger -Thiemer, Grodins, Jacquez, Segre and 
others, it is well established that there now exists sufficient physiological 
understanding to permit the construction of useful mathematical models 
of the distribution of specific drugs in the human body. In the standard 
fashion these models can be used to test and develop physiological concepts 
as well as to guide experimental procedures. Such mathematical models 
form the foundation for intensive investigations of cancer, viewed as a 
control, or rather, out-of-control, process, and other fundamental 
processes and, more generally, will provide the basis of one of the new 
fields of the future --mathematical physiology. 

An important application of pharmacokinetics is to the field of 
chemotherapy, in particular, to the avoidance of "side-effect" reactions 
of prescribed drugs such as antibiotics, kanamycin and others, or digoxin. 
The problem may be viewed as that of providing adequate concentrations of 
the drug at a particular site in the body without having undesirable amounts 
at another site. We have shown that quite simple compar tmental models 
can be used to good avail in the administration of digoxin, kanamyecin, and 
gentamycin; the joint work of J. Buell and Dr. R. Jelliffe of the USC 
Medical School should be referred to for detailed results. 

This activity in pharmacokinetics is being developed into a major 
program, working closely with USC Medical School faculty. Dr. R. Jelliffe 
and Dr. R. Maronde and his group. Both simple and complex mathematical 
and physiological models are being employed. 

Computer programs for suggested dosage regimens based upon drug 
kinetics are concurrently being used clinically through computer terminals 


- 83 - 



located at the Los Angeles County USC Medical Center, and other centers 
as well; see below. These programs include those for the cardiac 
glycosides; digitalis, digoxin, digitoxin- -including the metabolism of 
digitoxin, and ouabain. Concurrently, we are adding our latest computer 
models for use in antibiotic therapy for kanamycin and gentamycin which 
includes an inner ear compartment. Here the side effect is deafness. 

An integrated research program on therapy using these drugs under 
the supervision of Dr. R. Jelliffe is currently being undertaken at the 
LACH-USC medical facility, approved by the Research Council of the 
Department of Medicine. Volunteer patients are being used in a controlled 
hospital situation in order to validate computer-predicted serum level 
results with actual patient data. 

b. Nuclear Medicine 

In this area there are many important problems connected with both 
diagnosis and therapy, as well as many in the domain of operations 
research. One task is that of providing the physician with information 
concerning the presence and nature of a tumor, as well as of other 
abnormalities of an organ; a second is that of providing him with tech- 
niques for treating the tumorous condition. Both can involve the use of 
radiation. 

We know that early warning of cancer plays a vital role in successful 
treatment by means of surgery, radiology, and chemotherapy, all 
possibly combined in various ways. This means that vast screening pro- 
grams are required for minimally effective mass medical care. Valuable 
data can be obtained at the present using radioisotopes and scanning 
techniques. However, there are still numerous problems connected with 
the gathering and use of these data as well as with the efficient operation 
of giant programs of this nature. Questions of display of data and informa- 
tion become significant here. Statistical ideas and dynamic programming 
can be effectively used; the bios tatis -side of this work has been performed 
by the Friedland-Ueno-Vas ude van group. 

Let us briefly view some questions of therapy. In connection with the 
use of radiotherapy, it is essential to study the effects of radiation upon 
the human body. Once again, the side -effect problem is crucial. It is 
vital to ensure that the treatment does not harm the patient in the course 
of attempting to alleviate the original condition --the same problem is 
faced in chemotherapy. 

Since the human body consists of a set of inhomogeneous masses of 
irregular shape, the mathematical problems involved in ascertaining and 
controlling the effects of radiation are severe. Fortunately, our group 
possesses the advantage of a base of years of extensive research in 
astrophysics, transport theory, and nuclear physics, as well as the 
help of a number of distinguished experts in experimental and theoretical 
aspects of these areas; Professor S. Friedland, Professor R. Vasudevan, 
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and Professor S. Ueno, to cite a few. Invariant imbedding plays a major 
role here, as well as various other methods for solving integral equations 
and nonlinear equations. 

We are working closely with Dr. Siemsen at the USC Medical School, 
as well as with Dr. George's group In the Department of Radiology, LACH, 
and will actively pursue a number of other directions in radiotherapy. 

We have continued the work on the identification problem in cardiology, 
started some years ago with Drs. Collier and Selvester. Now that we have 
developed a powerful new mathematical technique in the theory of identifi- 
cation, differential quadrature, we intend to expand considerably the scope 
of the previous work built upon the basic physiological results of A. Scher. 

Research in neurophysiology involving cooperative effort of Dr. 

G. Moore of the Department of Electrical Engineering of USC, H. Sugiyama 
and S. Osaki has been continued and expanded. Again the newly developed 
technique of differential qudrature plays a significant role here. 

Work on large, complex control systems, with particular application 
to orthotics and prosthetics and stroke, has been continued and expanded. 
The development of parallel processing computers, (ILLIAC), may be the 
breakthrough required to make some existing ideas operational. 

c. Information Analysis 

This is a vast area which includes the major problem of making readily 
available to the physician the information he requires for effective decision- 
making. In order to make fruitful use of a limited supply of medical 
professionals, we wish to provide a variety of computer screening tech- 
niques based upon statistical and mathematical techniques, such as analysis 
of electrocardiograms and phonocardiograms, of tissue smears, and so on, 
as well as a number of new techniques for the display of data using computer 
graphics . 

This work is headed by S, Azen. Data are presented in one of many 
forms (e.g., histograms, scattergrams , etc.) on a CRT tube; subsamples 
are selected, manipulated, and analyzed. The researcher makes decisions 
along the way. Existing graphics and statistics programs (e. g. , the BMD 
package) will be utilized. Closely tied in with this is the task of developing 
a "biostatistics language": a language, simple to learn, that performs 
statistical manipulations. For example, "REGRESS Y ON X" performs a 
simple linear regression. 

A package of examples has been prepared which illustrates statistical 
concepts. This can be used as a valuable device to teach statistical tech- 
niques to medical people. A simple version is currently being used in 
teaching dental students at UCLA. In preparation is a somewhat more 
sophisticated package which explains sampling- -random or stratified, 
the likelihood ratio procedure, maximum likelihood estimation, and so 
forth. Also under consideration are the application, extension, and 
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comparison of graphic procedures based on non-parametric Fourier 
estimates of probability densities, cumulative distribution functions, 
survival curves , etc. 

Pattern recognition in general remains a stubborn field requiring a 
continuing effort. Despite its difficulty, it cannot be shunted aside since 
many kinds of pattern recognitions are absolutely essential for a health 
program for a nation of 200,000,000. A great deal can be done in 
connection with the identification of systems, as we shall indicate below. 

We expect that the parallel processing computer, such as ILLIAC-4 cited 
above, will allow dynamic programming to be used fruitfully in these areas. 

The modern theory of statistics, combined with computers, is essential 
for the storage and retrieval of data required for large-scale patient care. 
This is an essential part of the study of large health systems. 

All of this represents a sustained effort to analyze and improve the 
uses of computers in a hospital environment. It is part of a systematic 
program begun six years ago in conjunction with Dr. R. Maronde of LACH. 
We are also concerned with the training of the operators of computers and 
computing systems as well as the design of college curricula for these 
purposes . 

We have emphasized decision-making as an essential part of medical 
care. Decision-making involving quantitative effects, and thus numbers, 
can be handled by various traditional mathematical methods. However, 
medicine involves people, and thus human values. These cannot be treated 
wholly by numbers, a fact which renders decision-making in medical areas 
very much more difficult than engineering decision-making. In the 
operation of medical systems, however, we cannot avoid these more 
difficult problems, nor those connected with general man-machine systems. 
On-line decision-making plays a vital role here. We refer below to the 
use of simulation techniques we have developed to study these types of 
questions, making significant use of the computer in this activity. 

C. P. Kell is in charge of these areas, aided by consultants Dr. W. C. 

Hop good of UCLA, a practicing psychoanalyst, and Professor L. Zadeh 
of the Department of Engineering of the University of California at 
Berkeley, who is applying his new theory of fuzzy sets. There are many 
important applications to psychotherapy, both conceptually and operationally. 

e * Hospital and Community Clinic Administration 

In this area we apply operations analysis and system analysis tech- 
niques to the many problems arising in the areas of scheduling, 
communication, and organization. Many of these can be adequately 
treated by a combination of conventional mathematical theory and 
simulation. Dimensionality problems are severe, but again we expect 
that dynamic programming ideas will become much more effective with 
the new parallel-processing computers. S. Osaki has been working actively 
in these areas . 
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We are fortunate in being connected with Resthaven, a community 

mental health facility of moderate size, in the study of various problems 

involving the delivery of health services to the community, working 

closely with Dr. J. Krofcheck, Clincal Director, and his staff. 
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Objective: 

To study a wide range of basic and applied problems in biological dynamics, 
control and communication. The specific fields of interest are the respira- 
tory, cardiovascular, fluid -electrolyte, neural, and metabolic systems. 

Recent Progres s : 

Only a brief summary of major accomplishments is presented here. 
Details can be found in the references. 

a. Respiratory System 

We have studied problems in respirator y mechanics , cycle optimiza- 
tion, chemical control, and the hyperpnea of exercise. In the area of 
mechanics, we have explored the mechanisms limiting maximum expiratory 
flow rate. Most previous analyses of this problem have treated airway 
impedance in terms of "tubes" whereas in reality the tracheobronchial 
tree is largely a collection of "junctions 11 connected by segments so short 
as to hardly qualify as tubes. By a more rigorous application of aero- 
dynamic flow theory to a more realistic tracheobronchial flow net, we 
have shown that the phenomenon of "selective merging" at the junctions 
can account for 55% of the change in impedance during the expiratory flow 
maneuver, A probable cause of the remainder is junctional rather than 
tubular distortion, but this remains to be proven. 

We have explored the possibility that the mechanical respiratory 
system is regulated in a way which minimizes the work required to 
produce a given level of alveolar ventilation. This hypothesis was tested 
by analytically predicting optimal airflow patterns and comparing these 
to experimentally measured patterns available from the literature. Good 
agreement was found with many animal patterns and human patterns 
observed under stress conditions. Human airflow patterns at rest, 
however, could not be explained by this hypothesis. Instead, a criterion 
of minimum mean-squared volume acceleration was found to predict 
airflow patterns which match the resting human pattern. This criterion 
can be linked to gas transport efficiency. 

In tlie area of chemical control, we have explored the role of peripheral 


- 94 - 



and central receptors in the ventilatory responses to acute metabolic 
disturbances in acid-base balance. The responses of nembutalized dogs 
to infusions of 0. 5 MHC1 or 7% NaHCC >3 were studied with the carotid 
bodies intact or denervated and with Pa(C 02 ) c ^ am P e d constant through 

servo-control. It was concluded that both peripheral and central 
receptors contribute to the ventilatory responses under these conditions. 


That both receptors also contribute to the ventilatory response to 
sinusoidally varying in man has been supported by a computer 

simulation study. Using the mathematical model of Grodins, Buell, and 
Bart (1967), the responses to sinusoidal for several alternative 

controller hypotheses were explored, and then compared with experimental 
data in man taken from the literature (Bellville, 1969; Stoll, 1969). The 
best agreement was achieved when the controller equation included both a 
central and peripheral drive. This study also revealed that simulated 
responses to a sinusoidal input were much more sensitive to 

the form of the controller hypothesis than were responses to the more 
conventional step input. 


To help resolve the mystery of exercise hyperpnea, we have done 
some preliminary studies on the ventilatory responses to dynamic P a (co 2 ) 

P a (02) P at terns associated with the respiratory cycle. A new technique 


was developed to provide a positive inspired CO 2 load and a ''negative" 
inspired O 2 load which were independent of ventilation. That this procedure 
should produce P a (CO 2 ) " F a ( O 2 ) dynamic patterns similar to those of 

exercise was supported by a simulation study. Preliminary results 
in dogs indicated that the ventilatory response to this form of CO 2 -O 2 
loading resembles that to exercise. Recent results on work in this area 
can be found in References a, b, c and d. 


b. Cardiovascular System 

We have studied the baroceptor reflex in the fetal lamb, developed 
and applied a frequency response method for the measurement of 
peripheral vascular parameters, and developed a promising non-invasive 
method of estimating cardiac output from respiratory measurements. 

The baroceptor study was an extension of previous work carried out 
prior to the initiation of this program in which a mathematical model of 
the isolated, uncontrolled "cardiovascular plant" of the fetal lamb was 
developed. In the present study, the model was extended by adding a 
baroceptor control loop. The two major physiological questions asked 
in this study were (1) does a baroceptor reflex exist in the fetal lamb? 
and (2) if it does, what plant parameters does it operate on? Model 
simulations were compared with experimental data provided by Assali 
and Brinkman working at UCLA. It was concluded that a baroceptor 
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reflex does exist, and that it must operate in part upon the conductance of 
the placental circuit (References e and f). 

We have developed and applied a frequency response method for the 
measurement of mechanical and exchange parameters in the peripheral 
vascular bed of the dog. In this method, a small amplitude sinusoidal 
volume forcing is introduced at the arterial end of an isolated perfused 
hind-limb kept under isogravimetric conditions. The responses of 
arterial and venous pressures, which are also sinusoidal, are recorded 
and the procedure is repeated over a range of frequencies. Using these 
frequency response curves in conjunction with a mathematical model of 
the peripheral vascular bed, it is possible to obtain numerical values 
for seven parameters simultaneously. Thus arterial, venous, capillary, 
and tissue compliances, as well as pre- and post-capillary resistances 
were measured under a variety of conditions (carotid sinus stimulation, 
sympathetic stimulation, transfusion, epinephrine infusion) (Ref. g). A 
seventh parameter which has the dimensions of a resistance or conductance, 
may or may not be related to a capillary filtration coefficient. This point 
is under further study at present (References h and i). 

The parameter estimates mentioned above were based on a linear, 
lumped element model which assumed uniform tissue and capillary compartments. 
However, a CFC (capillary filtration coefficient) gradient has been 
reported from the arterial to the venous end of capillaries. Also, a 
hydrostatic gradient is known to be present along the capillary length. 

The impact of these characteristics on the estimate of a single CFC value 
was studied by computer simulation. It was found that gradients of CFC 
and pressure do not have significant effects. The estimate of a single 
CFC merely reflects the average CFC along the capillary length (Ref. j). 

While the frequency response method yields useful parameter 
estimates, its main disadvantage is the lengthy experimental period 
required to collect the necessary data. Up to one half hour of data are 
required during which stationarity of the preparation may be a problem. 

We have evaluated several alternative methods of estimation in an attempt 
to shorten the time required. Computer simulation results indicated that 
a modulating function approach was the most promising (Reference k). 

We have developed a non-invasive method for the measurement of 
cardiac output using only respirator y measurements and a mathematical 
model of the respiratory system. It is based upon the fact that cardiac 
output is a parameter of the respiratory system. A mathematical model 
of this system (Grodins and James, 1963) is simulated on a computer and 
forced (on-line) with a signal proportional to the subject's total ventilation. 

The output of the model is a signal representing end- tidal P C q . This 
computer generated Pet(CC> 2 ) is compared with the subject's ^ET(C0 2 ) 

as measured with a Beckman C0 2 analyzer. A criterion function is 
formed by squaring the difference between these two signals at the end 
of each breath and minimized by adjusting the model cardiac output using 
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a discrete s teepes t descent procedure. The method was tested in ten 
anesthetized dogs using the cardiogreen dye dilution technique as the 
primary standard of comparison. Model and dye estimates agreed 
within 10% for all experiments over a three-fold range of cardiac output. 
Over half of the estimates agreed within 4%. 

c. Fluid -Electrolyte Studies 

Studies in this area have been concerned with renal hemodialysis, 
dynamic interchanges of water and solutes between fluid compartments, 
and the renal regulation of bicarbonate balance. 

1. Renal hemodialysis. We have conducted investigations designed 
(1) to predict changes in body chemistry during artificial kidney (hemo- 
dialysis) treatments, (2) to optimize artificial kidney treatment in the 
sense of removing maximum amounts of chemical substances in a given 
time without producing gross imbalances in body chemical distributions, 
and (3) to quantitatively describe the dynamic and steady state rates of 
chemical removal by the artificial kidney machine. 

To aid the clinician in his prescription of artificial kidney treatment 
for an individual patient, we devised mathematical models which could 
accurately predict body fluid concentrations of urea and creatinine in a 
given patient knowing only his weight and the average blood flow into the 
kidney machine. The results were used to generate nomograms (1) from 
which the physician can rapidly make estimates of concentration changes 
for a given patient under treatment. 

We have done some theoretical studies on maximizing the rate of 
removal of urea from a patient while maintaining transcellular concentration 
gradients below specified levels. Presumably the maintenance of small 
concentration gradients would prevent the occurrence of the "Disequili- 
brium Syndrome" which comes about during rapid dialysis in patients. 

Those studied showed that the method of dynamic programming could be 
used to specify optimum removal rates under these constraints. 

In conjunction with the above studies, we have developed mathematical 
descriptions for the removal of chemical substances from the body by the 
artificial kidney (m). A transient analysis of these devices may prove 
important in studying the chemical changes in a patient during the initial 
minutes of hemodialysis where the blood flow is not constant and the 
chemical concentrations in the patient's blood are most rapidly changing. 

Dynamics of Com par tmental Exchange. We have constructed 
a mathematical model describing the dynamic interchanges of water and 
solute in an animal or human following rapid infusions of chemical 
solutions or other perturbations. The model uses the compar tmental 
concept and considers fluid interchanges between red cells, plasma, 
interstitial fluid and intracellular fluid to result from pressure and 
concentration differences between these fluids as well as from solvent- 
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solute coupling (see Reference n). 

Since the model includes the whole body rather than a single isolated 
perfused organ or tissue, the parameters are, for the most part, deter- 
mined indirectly by fitting to experimental data. It has thus been 
necessary to devise methods to measure a sufficient number of variables 
in the experimental animal to permit reliable estimates of these 
parameters. In particular, methods for measuring plasma volume 
dynamics were investigated (o). Two independent methods were used. 

The first depends upon the disappearance curve of radioactive albumin 
from the circulation, and the second upon dynamic measurements of red 
cell volume and hematocrit. Although each of the methods had been 
previously suggested in the literature, they had never been used simul- 
taneously or during rapid transients. 

Both methods gave very similar results for the transient responses 
of plasma volume after fluid infusion. These data together with transient 
changes of red cell volume and chemical concentrations have been used 
to establish whole body parameters for the model. We are currently 
evaluating the variance of these from animal to animal, as well as 
comparing them with estimates derived from isolated organ experiments. 

Models such as this can ultimately be used to assess the effects of 
fluid infusions and losses inhuman patients. 

3. Renal Control of Electrolyte Balance . We are currently 
conducting experiments to determine the most important factors respon- 
sible for the renal regulation of bicarbonate balance. In this current 
work we have centered our attention on the renal reabsorption of the 
bicarbonate ion, since it is through this process that the fine control of 
excretion is achieved. 

Previous workers have implicated a number of factors in the 
modification of bicarbonate reabsorption. These are: the plasma 
concentrations of hydrogen, bicarbonate and potassium ions; Pqq^J 

ratio of the reabsorbed to the filtered sodium; and the extracellular 
fluid volume (ECF). In all previous work, the changes in reabsorption due 
to a number of these factors were studied, and in every case the ECF was 
expanded in order to increase plasma bicarbonate concentration above the 
"renal threshold" as defined by Pitts and Lotspeich. Thus, the observed 
bicarbonate changes were always associated with simultaneous variations 
in a number of potential causative factors, and the individual effects of 
each could not be separated. Our experiments are unique in that we 
have raised plasma bicarbonate levels without expanding ECF. This 
was accomplished by exchange hemodialysis with an artificial kidney 
machine. The results indicate that ECF expansion is a primary but 
not the sole factor in the depression of bicarbonate reabsorption at high 
levels of plasma bicarbonate concentrations. 
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d. Neuronal Signals and Systems 

There are four principal areas of work that have been under active 
research since the beginning of the program project grant. 

1) Much work, especially in the early period of the grant, was 
spent completing a study of the basic theory and applications of a 
technique for the extracellular detection and characterization of synaptic 
interactions between neurons. The basic approach assumes that when 
two neurons are synaptically connected, the firing probability of the 
postsynaptic cell will be influenced by activity in the presynaptic 

cell, especially after short delays during which a postsynaptic potential 
appears. The technique depends on the calculation of the crosscorrelation 
histogram for the two cells, and we showed based on earlier studies 
(Moore, Perkel, Segundo, Ger stein, 1966, 1967) that the sepcific features 
of the crosscorrelogram reflected the existence as well as the amplitude, 
waveform, and polarity of the postsynaptic potential, and, in addition, 
reflected the statistical features of the firing pattern of both the pre- and 
post-synaptic cells and other neurons involved in the interaction. These 
techniques are now being used routinely in many studies conducted by 
other laboratories, as well as our own. 

2) In connection with these studies, we continued to apply general 
statistical computer techniques to the study of single and multiple spike 
train records. In one study, we used auto- and c ross correlation to 
characterize patterns of activity in single chemoreceptor fibers of the 
cat carotid body, and to show that these fibers exhibit a pronounced 
rhythmicity at the respiratory period and possibly a smaller cyclic 
pattern synchronous with the cardiac cycle* These techniques are 

now being extended to studies of the mechanism of transmitter release 
at the neuromuscular junction. Reference n relates to this work. 

3) As a result of these investigations, it became clear that some 
more powerful methods of analysis, modeling and simulation of neuronal 
systems would be desirable. We therefore developed a powerful compu- 
tational program for analysis and simulation of small neuronal networks 
so that the behavior of experimentally recorded cells and their simulated 
counterparts could be compared directly by statistical analysis. The 
unique feature of this system is that it was designed so that the parameters 
of the simulation program could be optimized with respect to the experi- 
mental data; this is not generally available and gives us enormous power 
in testing and validating proposed models of neuronal processes. 
Preliminary reports of this work have been made and will appear in 
extended form in several fortcoming papers. 

An additional development has been the mathematical solution of a 
general class of mathematical models for neuronal spike train activity. 
This work allows the determination of general constraints on the 
parameters of an observed neuron and will be the cornerstone of a 
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generalized computer program which can identify neuronal parameters 
from an observed interval histogram. 

4) Finally, we are presently conducting a number of studies whose 
aim is the development of more powerful and accurate models of single 
neuron processes. In particular, we are presently investigating models 
for the sodium pump, for transmitter synthesis and release, and have 
developed a general ionic flux model of the membrane, which, when 
combined with the Hodgkin-Huxley equations can account for an 
entire spectrum of passive sub- and s uprathr eshold benavior of a 
neuron (Reference q). 

e. Metabolic and Endocrine Systems 

This is the newst component of our research program , having only 
been underway since about the middle of 1970. Nevertheless, considerable 
progress has been made towards achieving the long term goals of this 
division. These goals include (1) the characterization of the organ system 
dynamics of liver, pancreas and adipose tissue in response to hormonal 
and substrate signals, (2) the definition of glucose and fatty acid (power) 
fluxes in the blood of unanesthetized animals as they adopt various 
metabolic states, (3) the identification of the corresponding informational 
fluxes associated with the power fluxes identified in (2) above. 

For the study of organ system dynamics, highly reliable, computer 
controlled perfusion systems have been developed and tested using the 
"pilot organ" approach. Since these are described in detail in the body 
of the proposal (pp. 172-176 for fat; pp. 179-182 for liver; pp. 191-193 
for pancreas), they need not be discussed further here. Some preliminary 
results using the pilot pancreas have recently been reported (r). In this 
study, temporal patterns of insulin secretion rate from the pilot pancreas 
in response to well-define variations in the pancreatic arterial concentra- 
tions of glucose and amino acids were determined. The results indicated 
that glucose and amino acids stimulate insulin secretion by distinct 
mechanisms and demonstrated the dependence of the amino acid 
stimulated response on arterial glucose concentration. 

In addition to these studies on isolated perfused organs or tissue, 
changes in both power fluxes (glucose and free fatty acids) and information 
fluxes (insulin) have been studied in conscious dogs. Some details of these 
studies are described on pp. 197-198 of the proposal, and some preliminary 
results have recently been reported (s). In the latter, temporal patterns 
of arterial plasma glucose concentrations during the resting state and 
during steady glucose infusions were studied in fasting conscious dogs 
over periods of thirty minutes to 1 1 hours. The results suggested that 
the glucose regulatory system may enter an oscillatory mode when 
stimulated by moderate infusions of glucose, but that this mode may not 
be detectable at rest or during low rates of infusion. 


- 100 - 



References and Recent Publications: 


a. Yamashiro, S. M. , and F. S. Grodins, "Optimal Regulation of 
R espiratory Airflow, " J*. Appl. Physiol. 30: 597-602, 1971. 

b. Yamashiro, S. M, , F. S. Grodins, and S. Sohrab, "Simulation of 
Rspiratory Responses to Sinusoidally Varying Inspired C02» " 1971 
San Diego Biomedical Symposium. 

c. Yamashiro, S. M. , F . S. Grodins, and M. Gordon, "Simulation of 
Alveolar I^COo Dynamics During External Metabolic CO 2 Loading, " 

24th ACEMB, 1971. 

d. Yamashiro, S. M. , F. S. Grodins, and M. Gordon, "Effect of 
Simulated Metabolic Loading on Respiration," Physiologist 14: 255, 1971. 

e. Bekey, G. A., W. Morrison, and N. S. Assali, "Mathematical 
Model of the Baroreceptor Reflex in the Full-Term Fetus, " Proc. 8th 
Int. Conf. on Medical and Biological Engineering, 20-24, 1969. 

f. Bekey, G. A., L. W. Morrison, C. R. Brinkman, and N. S. Assali, 
"Computer Simulation of Fetal Baroreceptor Function, " Computers 
and Biomedical Research, 3: 561-574, 1971. 

g. Sato, T., S. M. Yamashiro, and F. S. Grodins, "Measurement of 
Peripheral Vascular Properties by a Frequency Response Method, " 

Am. J. Physiol. 220: 1640-1650, 1971. 

h. Sato, T., S. M. Yamashiro, and F. S. Grodins, "Estimation of 
Capillary Hydrostatic Pressure and Filtration Coefficient Using 
Sinusoidal Perturbation of Isogra vimetric State," Fed. Proc. 30: 

661 Abs. , 1971. 

i. Yamashiro, S. M., T. Sato, and F. S. Grodins, "Gravimetric 
Determination of Tissue Exchange Parameters by Dynamic Analysis, " 

3rd Annual BMES Meeting, Abstract 17, 1972. 

j. Sato, T., S. M. Yamashiro, and F. S. Grodins, "Simulation of 
Local Control Dynamics of the Peripheral Vascular System, " 1971 
San Diego Biomedical Symposium. 

k. Yamashiro, S. M. , T. Sato, and F. S. Grodins, "Comparison of 
Parameter Identification Methods Applied to a Study of Peripheral 
Vascular Function, " Fed. Proc. 30: 698 Abs., 1971. 

l. Wolf, M. B. , B. H. Barbour, and P. D. Watson, "Nomograms for 

the Prediction of Changes in Patient Chemistries During Hemodialysis, l. 11 
Kidney Internat. 1: 126-130, 1972. 


- 101 - 


m. 


Buell, J. , R. Kalaba, and M. B. Wolf, "A New Approach to Transients 
in Countercurrent Flow Systems," Chem. Engr. Sci. (in press). 

n. Wolf, M. B. , "A Mathematical Model of Water and Solute Exchange, " 
Fed. Proc. 30: 703, 1971. 

o. Wolf, M. B. , "Plasma Volume Dynamics After Hypertonic Fluid 
Infusions in Nephrectomized Dogs, M Am. J. Physiol. 221: 13 92- 
1395, 1971. 

p. Gehrich, J. L., The Application of Statistical Methods to the 
Investigation of the Carotid Body Chemoreceptor Mechanism, 
University of Southern California, Ph.D. dissertation, June 1971. 

q. Kroin, J. , R. J. Sclabassi, and G. P. Moore, "Experimental Studies 
of Accommodation: A Theoretical Explanation," submitted for 
publication. 

r. Bergman, R. N. , "Glucose and Amino Acids : Contrasts and 

Interactions in Stimulating Insulin Secretion, " 3rd Annual BMES 
Meeting, Abstract 62, 1972. 

s. Marsh, D. J. , M. Ookhtens, R. N. Bergman, and F. E. Yates, 
"Fluctuations of the Glucose Concentration in. Arterial Plasma of 
Unanesthetized Dogs, With and Without Steady Infusions of Glucose, " 
3rd Annual BMES Meeting, Abstract 61, 1972, 


- 102 - 



